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The B-factories and Large Hadron Collider experiments have demonstrated the ability to observe 
and measure the properties of bottomonium mesons. In order to discover missing states it is useful 
to know their properties to develop a successful search strategy. To this end we calculate the masses 
and decay properties of excited bottomonium states. We use the relativized quark model to calculate 
the masses and wavefunctions and the ^Pq quark-pair creation model to calculate decay widths to 
open bottom. We also summarize results for radiative transitions, annihilation decays, hadronic 
transitions and production cross sections which are used to develop strategies to find these states. 
We find that the bb system has a rich spectroscopy that we expect to be substantially extended by 
the LHC and e''"e“ experiments in the near future. Some of the most promising possibilities at the 
LHC are observing the Xf)(i,2) (3P), Xi)(i.2) (dB) and ry,{SS) states in final states that proceed 

via radiative transitions through T(nS') intermediate states and l^Dj and 2^Dj into pT final 
states proceeding via GPj —^ and 2^Pj —> 2^51 intermediate states respectively. Some of the 

most interesting possibilities in e^e“ collisions are studying the l^Dj states via dy cascades starting 
with the T(3S) and the 3^Pj states in 'y'yp'^p~ final states starting with the T(4S') and proceeding 
via T{nS) intermediate states. Completing the bottomonium spectrum is an important validation 
of lattice QCD calculations and a test of our understanding of bottomonium states in the context 
of the quark model. 

PACS numbers: 12.39.Pn, 13.25.-k, 13.25.Gv, 14.40.Pq 


I. INTRODUCTION 

The Large Hadron Collider experiments have demon¬ 
strated that they can discover some of the missing bot¬ 
tomonium states. In fact, the first new particle discov¬ 
ered at the LHC was a 3P bottomonium state mm- 
The Belle H experiment at SuperKEKB will also offer 
the possibility of studying excited bottomonium states 
[3]. At the same time, lattice QCD calculations of bot¬ 
tomonium properties have advanced considerably in re¬ 
cent years HHi] so it is important to expand our exper¬ 
imental knowledge of bottomonium states to test these 
calculations. With this motivation, we calculate prop¬ 
erties of bottomonium mesons to suggest experimental 
strategies to observe missing states. The observation of 
these states is a crucial test of lattice QCD calculations 
and will also test the various models of hadron proper¬ 
ties. Some recent reviews of bottomonium spectroscopy 
are Ref. HIIQ]. 

We use the relativized quark model to calculate the 
masses and wavefunctions m- The mass predictions for 
this model are given in Section H. The wavefunctions are 
used to calculate radiative transitions between states, an¬ 
nihilation decays and as input for estimating hadronic 
transitions as described in Sections HI-V respectively. 
The strong decay widths to open bottom are described 
in Section VI and are calculated using the model 
[UlIIS] with simple harmonic oscillator (SHO) wavefunc¬ 
tions with the oscillator parameters, /3, found by fitting 
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the SHO wavefunction rms radii to those of the corre¬ 
sponding relativized quark model wavefunctions. This 
approach has proven to be a useful phenomenological tool 
for calculating hadron properties which has helped to un¬ 
derstand the observed spectra [I4UI8] . Additional details 
of the ^Pq model are given in the appendix, primarily so 
that the various conventions are written down explicitly 
so that the interested reader is able to reproduce our re¬ 
sults. 


We combine our results for the various decay modes 
to produce branching ratios (BR) for each of the bot¬ 
tomonium states we study. The purpose of this paper 
is to suggest strategies to find some of the missing bot¬ 
tomonium states in pp collisions at the LHC and in e“'"e“ 
collisions at SuperKEKB. The final missing input is an 
estimate of production rates for bottomonium states in pp 
and e+e” collisions. This is described in Section VH. We 
combine the cross sections with the expected integrated 
luminosities and various BR’s to estimate the number 
of events expected for the production of bottomonium 
states with decays to various final states. This is the 
main result of the paper, to identify which of the miss¬ 
ing bottomonium states are most likely to be observed 
and the most promising signal to find them. However, 
there are many experimental issues that could alter our 
conclusions so we hope that the interested reader can use 
the information in this paper as a starting point to study 
other potentially useful experimental signatures that we 
might have missed. In the final section we summarize the 
most promising signatures. 
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II. SPECTROSCOPY 


We calculate the bottomonium mass spectrum using 
the relativized quark model m- This model assumes a 
relativistic kinetic energy term and the potential incor¬ 
porates a Lorentz vector one-gluon-exchange interaction 
with a QCD motivated running coupling constant, as{r), 
and a Lorentz scalar linear confining interaction. The 
details of this model, including the parameters, can be 
found in Ref. [TT] (see also Ref. [TM22] L This is typi¬ 
cal of most such models which are based on some vari¬ 
ant of the Coulomb plus linear potential expected from 
QCD and often include some relativistic effects. The rel¬ 
ativized quark model has been reasonably successful in 
describing most known mesons and is a useful template 
against which to identify newly found states. However in 
recent years, starting with the discovery of the Dsj{2317) 
[53H5S] and Y(3872) states [2^, an increasing number of 
states have been observed that do not fit into this picture 
[27H5D] pointing to the need to include physics which has 
hitherto been neglected such as coupled channel effects 
m- As a consequence of neglecting coupled channel ef¬ 
fects and the crudeness of the relativization procedure 
we do not expect the mass predictions to be accurate to 
better than ~ 10 — 20 MeV. 

The bottomonium mass predictions for this model are 
shown in Fig. [T] These are also listed in Tables [I1|TT| 
along with known experimental masses and the effective 
SHO wavefunction parameters, /3. These, along with the 
masses and effective /3’s for the B meson states, listed in 
Table m are used in the calculations of the open bot¬ 
tom strong decay widths as described in Sec. lYll We note 
that the l^Pi and B meson states mix to form the 
physical IPi and \P[ states, as defined in Table III with 


a singlet-triplet mixing angle of Oip = —30.3° for bq or¬ 
dering. 

If available, the experimental masses are used as in¬ 
put in our calculations rather than the predicted masses. 
When the mass of only one meson in a multiplet has been 
measured, we shift our input masses for the remaining 
states using the measured mass and the predicted split¬ 
tings. Specifically, to obtain the r]b{n^So) masses (for 
n = 3,4,5,6) we subtracted the predicted n^Si — u^Sq 
splitting from the measured T{n^Si) mass [52]. For the 
Xt,(3P) states, we calculated the predicted mass differ¬ 
ences with respect to the Xf,i(3^Pi) state and subtracted 
them from the observed X&i(3^Pi) mass recently mea¬ 
sured by LHCb |55| . We used a similar procedure for the 
T(ID) mesons [52] as well as for the currently unobserved 
IP B mesons [52] listed in Table jlllj 


III. RADIATIVE TRANSITIONS 

Radiative transitions of excited bottomonium states 
are of interest for a number of reasons. First, they probe 
the internal structure of the states and provide a strong 
test of the predictions of the various models. Moreover, 


TABLE I: Masses and effective harmonic oscillator parameter 
values (/3) for S-, P- and D-wave Bottomonium mesons. 


Meson 

Mtheo (MeV) 

Me,p (MeV) 

P (GeV) 

T(D5i) 

9465 

9460.30 ± 0.26 “ 

1.157 

Vbil^So) 

9402 

9398.0 ± 3.2 “ 

1.269 

r{ 2 ^Si) 

10003 

10023.26 ±0.31 “ 

0.819 

rib{2^So) 

9976 

9999.0 ±3.5li;| “ 

0.854 

T(3®5i) 

10354 

10355.2 ±0.5 “ 

0.698 

ribiS^So) 

10336 

10337 

0.719 

r{PSi) 

10635 

10579.4 ± 1.2 “ 

0.638 

Vb{4^So) 

10623 

10567 ‘ 

0.654 

T(5®Si) 

10878 

10876 ± 11 “ 

0.600 

rib{5^So) 

10869 

10867 ‘ 

0.615 

T(6®5i) 

11102 

11019 ±8 “ 

0.578 

rib(6^So) 

11097 

11014 

0.593 

Xb2{l^P2) 

9897 

9912.21 ±0.26 ±0.31 “ 

0.858 

Xblil^Pl) 

9876 

9892.78 ±0.26 ±0.31 “ 

0.889 

X 6 o(l®Po) 

9847 

9859.44 ±0.42 ±0.31 “ 

0.932 

hbiPPl) 

9882 

9899.3 ± 1.0 “ 

0.880 

Xb2{2.^ P 2 ) 

10261 

10268.65 ±0.22 ±0.50 “ 

0.711 

Xbi{2^Pi) 

10246 

10255.46 ±0.22 ±0.50 “ 

0.725 

X6o(2^Po) 

10226 

10232.5 ±0.4 ±0.5 “ 

0.742 

hb{2^Pi) 

10250 

10259.8 ± 0.5 ± 1.1 “ 

0.721 

Xb2{3^P2) 

10550 

10528 ‘ 

0.640 

Xbl{S^Pl) 

10538 

10515.7t3:313:® ° 

0.649 

Xbo{3^Po) 

10522 

10500 

0.660 

hbiS^Pi) 

10541 

10519 * 

0.649 

X 62 ( 4 ^P 2 ) 

10798 

N/A 

0.598 

Xbi{4^Pi) 

10788 

N/A 

0.605 

Xbo{4^ Po) 

10775 

N/A 

0.613 

hb{4^Pi) 

10790 

N/A 

0.603 

Xb 2 {b^ P 2 ) 

11022 

N/A 

0.570 

Xbi{5^Pi) 

11014 

N/A 

0.576 

Xbo{5^Po) 

11004 

N/A 

0.585 

hb{5^Pi) 

11016 

N/A 

0.575 

Tsil^Da) 

10155 

10172 ‘ 

0.752 

T2(1^D2) 

10147 

10163.7 ± 1.4 “ 

0.763 

Ti(DDi) 

10138 

10155 * 

0.776 

r/ 62 (l^D 2 ) 

10148 

10165 ‘ 

0.761 

T3(23D3) 

10455 

N/A 

0.660 

r 2 { 2 ^D 2 ) 

10449 

N/A 

0.666 

Ti(23Di) 

10441 

N/A 

0.672 

r]b2{2} D 2 ) 

10450 

N/A 

0.665 


10711 

N/A 

0.609 

r2{3^D2) 

10705 

N/A 

0.613 

Ti(33Di) 

10698 

N/A 

0.618 

7762(3^112) 

10706 

N/A 

0.612 

T 3 ( 434 ) 3 ) 

10939 

N/A 

0.577 

r 2 ( 4 ^ 02 ) 

10934 

N/A 

0.580 

Ti(43Di) 

10928 

N/A 

0.583 

7762(4'^ D2) 

10935 

N/A 

0.579 


“Measured mass from Particle Data Group m- 
^Using predicted multiplet mass splittings with measured mass as 
described in Sec.ITT| 

“Measured mass^om LHCb m- 
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FIG. 1: The bb mass spectrum as predicted by the relativized quark model m- 


for the purposes of this paper they provide a means of 
accessing bb states with different quantum numbers. Ob¬ 
servation of the photons emitted in radiative transitions 
between different bb states was in fact how the 3P bb 
state was observed by the ATLAS collaboration mm 
and subsequently by LHCb [321 IMj- El radiative partial 
widths of bottomonium are typically 0(1 — 10) keV so 
can represent a significant BR for bb states that are rela¬ 
tively narrow. As we will see, a large number of bb states 
fall into this category. With the high statistics available 
at the LHC it should be possible to observe some of the 
missing bb states with a well constrained search strategy. 
Likewise, SuperKEKB can provide large event samples of 
the T(3S') and T(45') and possibly the T(55') and T(6S') 
which could be used to identify radially excited P and 
D-wave and other high L states. e“''e“ collisions at Su¬ 
perKEKB could also produce the T(l^Zli) and T(2^Zli) 
directly which could be observed by Belle II in decay 
chains involving radiative transitions. 

We calculate the El radiative partial widths using [35] 

Tin^^+^Lj ^ ( 1 ) 


= 3 ^ Cfi5L.L'±i\ ( 7 /’/ |r I V> 2 ) P 

where the angular momentum matrix element is given by 

Cp = max(L, L'){2J' + 1) | ^ | ^ (2) 

and {)))} is a 6-j symbol, Cf, = —1/3 is the 6-quark charge 
in units of |e|, a is the fine-structure constant, is the 
photon energy and is the transition matrix el¬ 

ement from the initial state ipi to the final state ipj. For 
these initial and final states, we use the relativized quark 
model wavefunctions m- The El radiative widths are 
given in Tables jTVj-jXXTTTj along with the matrix elements 
so that the interested reader can reproduce our results. 
The initial and final state masses are also listed in these 
tables where Particle Data Group (PDG) [33] masses are 
used when the masses are known. For unobserved states 
the masses are taken from the predicted values in Ta¬ 
bles |lj|ll] except when a member of a multiplet has been 
observed. In this latter case the mass used was obtained 
using the procedure described in Section |Tl| 
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TABLE II: Masses and effective harmonic oscillator parame¬ 
ter values (/3) for F- and G-wave Bottomonium mesons. 


Meson 

Mth.o (MeV) (MeV) 5 (GeV) 

Xi>4(l'P4) 

10358 

N/A 

0.693 

Xi>3(l'P3) 

10355 

N/A 

0.698 

Xb2{l^F2) 

10350 

N/A 

0.704 

hbsil^Fs) 

10355 

N/A 

0.698 

Xb4{2^F4) 

10622 

N/A 

0.626 

Xi,3(2^P3) 

10619 

N/A 

0.630 

Xb2{2^F2) 

10615 

N/A 

0.633 

hb3{2^F3) 

10619 

N/A 

0.629 

Xi,4(3^P4) 

10856 

N/A 

0.587 

XbsiS^Fa) 

10853 

N/A 

0.590 

Xm(3^P2) 

10850 

N/A 

0.592 

hbsiS^Fs) 

10853 

N/A 

0.589 

T5(1"G5) 

10532 

N/A 

0.653 

T4(1®G4) 

10531 

N/A 

0.656 

T3(1"G3) 

10529 

N/A 

0.660 

T}b4:{^^ G 4 ) 

10530 

N/A 

0.656 

T5(23G5) 

10772 

N/A 

0.602 

T4(23G4) 

10770 

N/A 

0.604 

T3(23G3) 

10769 

N/A 

0.606 

7754 (2^ G 4 ) 

10770 

N/A 

0.604 


An interesting observation is that the El transitions 
35' —>■ IP are highly suppressed relative to other El 
transitions [3H] (see Ref. m for a detailed discussion). 
Grant and Rosner |37! showed this to be a general prop¬ 
erty of El transitions, that El transitions between states 
that differ by 2 radial nodes are highly suppressed rela¬ 
tive to the dominant El transitions and are in fact zero 
for the 3-dimensional harmonic oscillator. As a conse¬ 
quence, these radiative transitions are particularly sen¬ 
sitive to relativistic corrections m- We found that this 
pattern was also apparent for similar transitions of the 
type In, Z) —t In — 2, Z ± 1) such as 55 —)■ 3P, 3P —)■ ID, 
4P ^ 25, 4D 2P, 3P> ^ IP etc. 

Ml transition rates are typically weaker than El rates. 
Nevertheless they have been useful in observing spin- 
singlet states that are difficult to observe in other ways 
[3H1I3S]- The Ml radiative partial widths are evaluated 
using BO] 


Tin^^+^Lj 


,2S' + 1 j 
n Lj, 


■l) 


(3) 


4aem 2J' + 1 
3ml 2L + 1 


Ss,S'±i\{i^f\jo{kr/2)\'tlji)\^ 


where jo (cc) is the spherical Bessel function and the other 
factors have been defined above. As with the El transi¬ 
tions, we use the relativized quark model wavefunctions 
m for the initial and final states. 

The partial widths and branching ratios for the Ml ra¬ 
diative transitions are listed in Tables[rVl|XXlII as appro¬ 


priate. For comparison, other calculations of bb radiative 
transitions can be found in Ref. [ssmiHis] 


IV. ANNIHILATION DECAYS 

Annihilation decays into gluons and light quarks make 
significant contributions to the total widths of some bb 
resonances. In addition, annihilation decays into lep¬ 
tons or photons can be useful for the production and 
identification of some bottomonium states. For exam¬ 
ple, the vector mesons are produced in e+e” collisions 
through their couplings to e“'"e“. Annihilation decay 
rates have been studied extensively using perturbative 
QCD (pQCD) methods [401HBH55] . The relevant formu¬ 
las for 5- and P-wave states including first-order QCD 
corrections (when they are known) are summarized in 
Ref. [S2] . Expressions for D- and P-wave decays are given 
in Refs. [SSJ [SB] and Refs. [SH |S7] respectively. The ex¬ 
pression for ^Di —>• e+e“ including the QCD correction 
comes from Ref. jSS] . Ackleh, Barnes and Close [S3| give 
a general expression for singlet decays to two gluons. A 
general property of annihilation decays is that the decay 
amplitude for a state with orbital angular momentum 
Z goes like where is the Z-th derivative 

of the radial wavefunction. P’1*) is typically 0{1) so for 
bottom quark masses the magnitude of the annihilation 
decay widths decreases rapidly as the orbital angular mo¬ 
mentum of the bottomonium state increases. Expressions 
for the decay widths including first-order QCD correc¬ 
tions when known are summarized in Table IXXIVI To 
obtain our numerical results for these partial widths we 
take the number of light quarks to be n/ =4, assumed 
ruh = 4.977 GeV, as ~ 0.18 (with some weak mass de¬ 
pendence), and used the wavefunctions found using the 
model of Ref. [n] as described in Section [TI] 

Considerable uncertainties arise in these expressions 
from the model-dependence of the wavefunctions and 
possible relativistic and QCD radiative corrections (see 
for example the discussion in Ref. m)- One example is 
that the logarithm evident in some of these formulas is 
evaluated at a rather arbitrarily chosen scale, and that 
the pQCD radiative corrections to these processes are of¬ 
ten found to be large, but are prescription dependent and 
so are numerically unreliable. As a consequence, these 
formulas should be regarded as estimates of the partial 
widths for these annihilation processes rather than pre¬ 
cise predictions. The numerical results for partial widths 
for the annihilation processes are included in Tables jlVf 

ixxni 


V. HADRONIC TRANSITIONS 

Hadronic transitions between quarkonium levels are 
needed to estimate branching ratios and potentially 
offer useful signatures for some missing bottomonium 
states. There have been numerous theoretical estimates 
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TABLE III: Masses and effective /3 values for B mesons used in the calculations of bottomonium strong decay widths. The 
physical IPi and IPi states are mixtures of l^Pi and l^Pi with singlet-triplet mixing angle 9ip = —30.3° for bq ordering. 
Where two values of f3 are listed, the first (second) value is for the singlet (triplet) state. 


Meson 

State 

Mtheo (MeV) (MeV) 

P (GeV) 


I'S'o 

5312 

5279.26 ±0.17 ° 

0.580 

B° 

l'5o 

5312 

5279.58 ±0.17 “ 

0.580 

B* 

l^Si 

5371 

5325.2 ±0.4 “ 

0.542 

B(l"Po) 

l^Po 

5756 

5702 

0.536 

B(1A) 

cos0ip{P Pi) + sin0ip(l®Pi) 

5777 

5723.5 ±2.0 “ 

0.499, 0.511 

B{lPl) 

— sin6ip(l^Pi) -1- cos0ip(l^Pi) 

5784 

5730 '' 

0.499, 0.511 

B{I^P2) 

1®P2 

5797 

5743 ± 5 “ 

0.472 

Bs 

l'5o 

5394 

5366.77 ±0.24 “ 

0.636 

Bt 

I’^S'i 

5450 

5415.4^^;^ “ 

0.595 


“Measured mass from Particle Data Group Ea- „ 

*’Input mass from predicted mass splittings, as described in Sec.[n] 


of hadronic transitions over the years [5^72) . In some 
cases the estimates disagree by orders of magnitude [53] . 
Hadronic transitions are typically described as a two-step 
process in which the gluons are first emitted from the 
heavy quarks and then recombine into light quarks. A 
multipole expansion of the colour gauge held is employed 
to describe the emission process where the intermediate 
colour octet quarkonium state is typically modeled by 
some sort of quarkonium hybrid wavefunction [601 172] . 
An uncertainty in predictions arises from how the re- 
hadronization step is estimated. To some extent this 
latter uncertainty can be reduced by employing the mul¬ 
tipole expansion of the colour gauge helds developed by 
Yan and collaborators [551 - 153] together with the Wigner- 
Eckart theorem to estimate the El-El transition rates 

m- 

In addition to EI-EI transitions such as —)■ 

there will be other transitions such as —)■ -|- rj, 

which goes via Ml-Ml & E1-M2 multipoles and spin-hip 
transitions such as —)■ which goes via El-Ml 

[60| . These transitions are suppressed by inverse powers 
of the quark masses and are expected to be small com¬ 
pared to the El-El and electromagnetic transitions. As 
a consequence, we will neglect them in our estimates of 
branching ratios. We note however, that in certain sit¬ 
uations they have provided a pathway to otherwise dif- 
hcult to observe states such as the he and /if, [TS] [74] 
and have played an important role in these states’ dis¬ 
coveries [75] [75]. Another example of a higher multipole 
transition is Xbi,b2{‘2P) —f a;T(lS') [77] which proceeds 
via three El gluons although it turns out that this par¬ 
ticular example has a larger branching ratio than the 
2P — )■ IP -I- TTTT transition [53] . 

The differential rate for El-El transitions from an ini¬ 
tial quarkonium state to the final quarkonium state 
$, and a system of light hadrons, /i, is given by the ex¬ 
pression [551155] : 


dr 

dA42 




$ -|- /i] — {2J P 1) ^ ^ 




Afc(/, i) 


( 4 ) 


where P, £ are the orbital angular momentum and P, J 
are the total angular momentum of the initial and hnal 
states respectively, s is the spin of the QQ pair, A4^ is 
the invariant mass squared of the light hadron system, 
and Ak{£', t) are the reduced matrix elements. For the 
convenience of the reader we give the expressions for the 
transition rates in terms of the reduced matrix elements 
in Table XXV The magnitudes of the Afc(t", t) are model 
dependent with a large variation in their estimates. The 
Afe(P, €) are a product of a phase space factor, overlap 
integrals with the intermediate hybrid wavefunction and 
a fitted constant. There is a large variation in the pre¬ 
dicted reduced rates. For example, for the transition 
l^P>i —1^5'i -I- TTTT, estimates for ^2(2,0) differ by al¬ 
most three orders of magnitude [55] [53 [751171] • In an 

attempt to minimize the theoretical uncertainty we esti¬ 
mate the reduced matrix elements by rescaling measured 
transition rates by phase space factors and interquark 
separation expectation values. While imperfect, we hope 
that this approach captures the essential features of the 
reduced matrix elements and gives a reasonable order of 
magnitude estimate of the partial widths. In the soft- 
pion limit the Ai contributions are suppressed so, as is 
the usual practice, we will take Ai{t', £) = 0 [55] (see 
also Ref. [75]) so that in practice only Aq{£', £) and/or 
A 2 {£', £) will contribute to a given transition. The Aq 
and A 2 amplitudes have phase space integrals of the form 

[53: 


G= -^TT^ 
4 M, . 


dM^^K I 1 - 


4toJ 

M2 

■‘-TTTT 



2 \2 


-2mi) 


(5) 


k £' £ 
s J J' 


2 


and 
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TABLE IV: Partial widths and branching ratios for strong and electromagnetic decays and transitions for the IS' and 2S 
bottomonium states. The state’s mass is given in GeV and is listed below the state’s name in column 1. Column 4, labelled M 
gives the matrix element appropriate to the particular decay; for S-wave annihilation decays A4 designates 4'(0) = i?(0)/\/4w 
representing the wavefunction at the origin and for radiative transitions the El or Ml matrix elements are ("i/’/klV’i) (GeV“^) 


i{kr/2)\ipi) respectively. Details of the calculations are given in the text. 

Initial 

Final 

Mf 

M 

Predicted 

Measured 

state 

state 

(GeV) 


Width (keV) 

BR (%) 

Width (keV) 

BR (%) 

T(135'i) 

l+£- 


0.793 

1.44 

2.71 

1.34 ± 0.04 

2.48 ±0.05“ 

9.460 “ 

ggg 


0.793 

47.6 

89.6 

44.1 ± 1.1 

81.7 ±0.7“ 


igg 


0.793 

1.2 

2.3 

1.2 ±0.3 

2.2 ±0.6“ 


777 


0.793 

1.7 X 10"® 

3.2 X 10"® 




ryi,(l^S'o)7 

9.398'" 

0.9947 

0.010 

0.019 




Total 



53.1 

100 

54.02 ± 1.25 ■" 


Vbi^^So) 

gg 


1.081 

16.6 MeV 

~ 100 



9.398 “ 

77 


1.081 

0.94 

0.0057 




Total 



16.6 MeV 

100 

10.8t^;^ MeV "" 


T(235'i) 

£+£- 


0.597 

0.73 

1.8 

0.62 ±0.06 

1.93 ±0.17 

10.023 “ 

ggg 


0.597 

26.3 

65.4 

18.8 ± 1.6 

58.8 ± 1.2 


igg 


0.597 

0.68 

1.7 

2.81 ± 0.42 

8.8 ± 1.1 


777 


0.597 

9.8 X 10"® 

2.4 X 10"® 




X62(1^E2)7 

9.912'" 

-1.524 

1.88 

4.67 

2.29 ± 0.22 

7.15 ±0.35“ 


Xf,i(l^Pi)7 

9.893'" 

-1.440 

1.63 

4.05 

2.21 ± 0.22 

6.9 ±0.4“ 


X6o(1^Eo)7 

9.859'" 

-1.330 

0.91 

2.3 

1.22 ±0.16 

3.8 ±0.4“ 


77i,(2^S'o)7 

9.999'" 

0.9924 

5.9 X 10"4 

1.5 X 10"® 




i7i)(l^*S'o)7 

9.398'" 

0.0913 

0.081 

0.20 

0.012 ±0.004 

(3.9 ± 1.5) X 10"® “ 


T(l35'i)7r7r 



8.46"" 

21.0 

8.46 ±0.71 

26.45 ± 0.48“ 


Total 



40.2 

100 

31.98 ±2.63 


VbC2^So) 

gg 


0.718 

7.2 MeV 

~ 100"" 



9.999 “ 

77 


0.718 

0.41 

5.7 X 10"® 




hi,(l4Pi)7 

9.899'" 

-1.526 

2.48 

0.034 




T(i3S'i)7 

9.460 

-0.0610 

0.068 

9.4 X 10"4 




J7i,(l^S'o)7r7r 



12.4 

0.17 




Total 



7.2 MeV 

100 

< 24 MeV"" 



“PDG Ref. [32]. 


H = 


1 M 


20 Mi 




dM^^K 


1 - 


4m^ 

M2 


1/2 


1 + 


2 


+ 


15M4 




( 6 ) 


respectively where 


tween decay rates for two bb states is given by 


(7) 

The amplitudes for El-El transitions depend quadrati- 
cally on the interquark separation so the scaling law be- 


r($i) (r2($i))" 

r($2) (r2($2))2- 


Because each set of transitions uses different experimen¬ 
tal input we will give details of how we rescale the 
sector by sector in the following subsections and 
give the predicted partial widths in the summary tables. 
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TABLE V: Partial widths and branching ratios for strong and electromagnetic decays and transitions for the 3S states. See 
the caption to Table [IV| for further details. 


Initial 

state 

Final 

state 

Mj 

(GeV) 

M 

Predicted 

Width (keV) BR (%) 

Measured 

Width (keV) BR (%) 

T(3^Si) 

£+£- 


0.523 

0.53 

1.8 


0.44 ± 0.06 

2.18 ±0.21 

10.355 “ 

999 


0.523 

19.8 

67.9 


7.25 ±0.84 

35.7 ±2.6 


799 


0.523 

0.52 

1.8 


0.20 ±0.04 

0.97 ±0.18 


777 


0.523 

7.6 X 10"® 

2.6 X 10” 

-5 




Xb2{2^P2)7 

10.269“ 

-2.446 

2.30 

7.90 


2.66 ±0.40 

13.1 ± 1.6“ 


Xbi{2^Pi)7 

10.255“ 

-2.326 

1.91 

6.56 


2.56 ±0.34 

12.6 ± 1.2“ 


Xbo{2^Po)7 

10.232“ 

-2.169 

1.03 

3.54 


1.20 ±0.16 

5.9 ±0.6“ 


Xf>2(l®P2)7 

9.912“ 

0.096 

0.45 

1.5 


0.20 ±0.03 

0.99 ±0.13“ 


Xi>i(l®-Pi)7 

9.893“ 

0.039 

0.05 

0.2 


0.018 ±0.010 

0.09 ±0.05“ 


Xi>o(l®Po)7 

9.859“ 

-0.028 

0.01 

0.03 


0.055 ±0.010 

0.27 ±0.04“ 


Y7i,(3^So)7 

10.337“ 

0.9920 

2.5 X lO"'* 

8.6 X 10” 

-4 




7?i,(2^So)7 

9.999“ 

0.1003 

0.19 

0.65 


< 0.12 

< 0.062 at 90% C.L.“ 


7?i,(RSo)7 

9.398“ 

0.0427 

0.060 

0.20 


0.01 ±0.002 

0.051 ±0.007“ 


T(l'^S'i)YrYr 



1.34“ 

4.60 


1.335 ±0.125 

6.57 ±0.15“ 


T(235'i)YrYr 



0.95“ 

3.3 


0.949 ± 0.098 

4.67 ±0.23 “ 


Total 



29.1 

100 


20.32 ± 1.85 


VbiS^So) 

99 


0.601 

4.9 MeV 

~ 100 




10 . 337 '' 

77 


0.601 

0.29 

5.9 X 10” 

-3 





10.260“ 

-2.461 

2.96 

0.060 





/ii,(l'Pi)7 

9.899“ 

0.1235 

1.3 

0.026 





T(235'i)7 

10.023“ 

-0.0484 

9.1 X 10"® 

1.8 X 10” 

-4 




T(1'^5'i)7 

9.460“ 

-0.031 

0.074 

1.5 X 10” 

-3 




?7i,(l^So)YrYr 



1.70 ±0.12 

3.5 X 10' 

-2 




77i,(2^So)YrYr 



1.16 ±0.10 

2.4 X 10” 

-2 




Total 



4.9 MeV 

100 





“PDG Ref. [32]. 

^Using predicted 3^Si — 3^So splitting and measured 3®Si mass. 


A. n'^So —t n^So +-TTYr 


The n'^So —>■ ri^So + tttt partial widths are found by 
rescaling the measured n'^Si —>■ n^Si +7r7r partial widths 
[32]. The n'S —>■ nS + tttt transitions are described by 
Ao(0, 0) amplitudes so that r(n'^S'o —?> n^5'o + tttt) are 
given by: 


r(n''5'o) = 


(r2(n'i5o))2 G{n 

■ X 


/I 


S'o —t n^Somr) 


r‘^{n'^Si))‘^ G(n'^Si n^Siinr) 


xr(n'^S'i). 

(9) 

The hadronic transition partial widths for the n'^So 
states are given in Tables Iv]|vl for n' = 2, 3,4. 

We do not make predictions for the rj^ibS) state as 
the measured hadronic transition rates for the T(5iS') are 
anomalously large and inconsistent with other transitions 
between 5'-waves Haj. This has resulted in speculation 
that the T{5S) is mixed with a hybrid state leading to its 
anomalously large hadronic transition rates |80j , contains 


a sizable tetraquark component lEiiiHa or is the conse¬ 
quence of rescattering [53]. See also Ref. [31 [Hj- 

It could also be the result of a large overlap with the in¬ 
termediate states. This subject needs a separate more de¬ 
tailed study which lies outside the present work. We also 
do not include hadronic transitions for the 65" states as 
there are no measurements of hadronic transitions orig¬ 
inating from the 6^Si state and in any case, the total 
widths for the 6S states are quite large so that the BR’s 
for hadronic transitions would be rather small. 


B. n'Pj —>■ nPj + TTYT 


All n'Pj —>■ nPj + TTTT transitions can be expressed 
in terms of Ao(l,l) and ^ 2 ( 1 ,!) where we have taken 
Ai(l, 1) = 0. The expressions relating the various partial 
widths in terms of these reduced amplitudes are summa¬ 


rized in Table XXV We can obtain Ao(l, 1) and ^ 2 ( 1 ,1) 













TABLE VI: Partial widths and branching ratios for strong and electromagnetic decays and transitions and OZI allowed strong 
decays for the dS and 5S bottomonium states. Details of the OZI allowed decay amplitudes are described in the appendix. See 
the caption to Table [IV] for further details. 


Initial 

Final 

Mf 

M 

Predicted 

Measured 

state 

state 

(GeV) 


Width (keV) 

BR (%) 

Width (keV) 

BR (%) 

T{4^Si) 

e+£- 


0.459 

0.39 

1.8 X 10"*^ 

0.32 ±0.04“ 

(1.57 ±0.08) X 10"®“ 

10.579“ 

ggg 


0.459 

15.1 

0.0686 




-ygg 


0.459 

0.40 

1.8 X 10"*^ 




777 


0.459 

6.0 X 10“® 

2.7 X 10"® 




Xb2{S^P2h 

10.528'* 

-3.223 

0.82 

3.7 X 10"® 




Xbi{3^Pih 

10.516"* 

-3.072 

0.84 

3.8 X 10"® 




Xbo{3^Poh 

10.500"* 

-2.869 

0.48 

2.2 X 10"® 




r{l^Si)n+TV- 



1.66*’ 

7.54 X 10"® 

1.66 ±0.24*’ 

(8.1 ± 0.6) X 10"®“ 


T{2^Si)n+Tv- 



1.76*’ 

8.00 X 10"® 

1.76 ±0.34*’ 

(8.6 ± 1.3) X 10"®“ 


BB 



22.0 MeV 

~ 100 


> 96 at 95% C.L. “ 


Total 



22.0 MeV 

100 

20.5 ± 2.5 MeV“ 


»7b(4'So) 

gg 


0.500 

3.4 MeV 

~ 100 



10.567“ 

77 


0.500 

0.20 

5.9 X 10"® 




hbiS^Pih 

10.519 "* 

-3.238 

1.24 

3.6 X 10"® 




?76(l^S'o)7r+7r“ 



2.03 ±0.29 

6.0 X 10"® 




?76(2^5o)7r+7r“ 



1.90 ±0.36 

5.6 X 10"® 




Total 



3.4 MeV 

100 



T(5®S'i) 

£+£- 


0.432 

0.33 

1.2 X 10"® 

0.31 ±0.23“ 

(5.6 ±3.1) X 10""*“ 

10.876“ 

ggg 


0.432 

13.1 

4.78 X 10"® 




Xb2{4^ P2)'y 

10.798 

-3.908 

4.3 

1.6 X 10"® 




Xm(4®-Pi)7 

10.788 

-3.724 

3.4 

1.2 X 10"® 




Xbo{^^Po)x 

10.775 

-3.483 

1.5 

5.5 X 10"® 




Xi)2(3®P2)7 

10.528 "* 

0.131 

0.42 

1.5 X 10"® 




Xm(3®Pi)7 

10.516 "* 

0.0020 

6.2 X 10“® 

2.3 X lO"’’ 




Xi>o(3®Po)7 

10.500 "* 

-0.156 

0.15 

5.5 X 10""* 




BB 



5.35 MeV 

19.5 

3.0 ± 1.6 MeV*" 

5.5 ± 1.0“ 


BB* 



16.6 MeV 

60.6 

7.5 ± 3.9 MeV*” 

13.7 ± 1.6“ 


B*B* 



2.42 MeV 

8.83 

21± 11 MeV*" 

38.1 ±3.4“ 


BsBs 



0.157 MeV 

0.573 

0.3 ±0.3 MeV*" 

0.5 ±0.5“ 


BsB* 



0.833 MeV 

3.04 

0.74 ±0.42 MeV*" 

1.35 ±0.32“ 


D* D* 

^ s ^ s 



2.00 MeV 

7.30 

9.7 ± 5.1 MeV*" 

17.6 ±2.7“ 


Total 



27.4 MeV 

100 

55 ± 28 MeV“ 


Vb{5^So) 

gg 


0.464 

2.9 MeV 

13 



10.867“ 

77 


0.464 

0.17 

7.4 X 10""* 




hi,(4iPi)7 

10.790 

3.925 

7.5 

3.3 X 10"® 




hi,(3iPi)7 

10.519 "* 

0.162 

1.1 

4.8 X 10"® 




BB* 



13.1 MeV 

57.0 




B*B* 



0.914 MeV 

3.97 




BsBl 



0.559 MeV 

2.43 




D* D* 

^ S ^ S 



5.49 MeV 

23.9 




Total 



23.0 MeV 

100 




as appropriate. 

“From PDG Ref. 1321 . 

^Found by combining the PDG BR with the PDG total widths 
for the T(43 Si) or r{5^Si) [32] 

“Using predicted n^Si —n^So splitting and measured n^Si mass. 

‘*3^Pi from LHGb 1331 and 3^P2, 3^Po and 3^Pi using predicted 
splittings with 3®Pi measured mass. 











9 


TABLE VII: Partial widths and branching ratios for strong and electromagnetic decays and transitions and OZI allowed strong 
decays for the 6S states. Details of the OZI allowed decay amplitudes are described in the appendix. See the caption to 
Table HVl for further details. 


Initial 

Final 

Mf 

M 

Predicted 

Measured 

state 

state 



Width (keV) 

BR (%) 

Width (keV) 

BR (%) 

T{6^Si) 

£+r 


0.396 

0.27 

8.0 X 10“'‘ 

0.13 ±0.05“ 

(1.6 ±0.5) X 10-'‘“ 

11.019“ 

ggg 


0.396 

11.0 

0.0324 




X62(5^P2)7 

11.022 


below threshold 

- 




X6i(5^Pi)7 

11.014 

-4.282 

8.3 X 10"^ 

2.4 X 10“® 




X6o(5®Po)7 

11.004 

-3.992 

6.4 X 10"® 

1.9 X 10“® 




X62(4®P2)7 

10.798 

0.116 

8.5 X 10"^ 

2.5 X 10“'‘ 




Xbi{4^Pih 

10.788 

-0.054 

0.012 

3.5 X 10“® 




X6o(4®Po)7 

10.775 

-0.244 

0.1 

3 X 10-^ 




X62(3®P2)7 

10.528 “ 

0.089 

0.53 

1.6 X 10“® 




X6i(3®Pi)7 

10.516 “ 

0.100 

0.43 

1.3 X 10“® 




X6o(3®Po)7 

10.500 “ 

0.115 

0.21 

6.2 X 10-* 




BB 



1.32 MeV 

3.89 




BB* 



7.59 MeV 

22.4 




BP(lPi) 



7.81 MeV 

23.0 




BP(lPi') 



10.8 MeV 

31.8 




B*B* 



5.89 MeV 

17.4 




Bs Bs 



1.31 

3.86 X 10“^ 




Bs Bg 



0.136 MeV 

0.401 




D* D* 

-^s 



0.310 MeV 

0.914 




Total 



33.9 MeV 

100 

79 ± 16 MeV“ 


76(6'So) 

gg 


0.410 

2.2 MeV 

16 



11.014* 

77 


0.410 

0.14 

1.0 X 10“® 




hbi5^Pi)x 

11.016 


below threshold 

- 




hbii^Pih 

10.790 

0.136 

0.22 

1.6 X 10“® 




hb{3^Pih 

10.519“ 

0.123 

1.8 

1.3 X 10-2 




BB* 



8.98 MeV 

66.0 




BB{l^Po) 



0.745 MeV 

5.48 




B*B* 



1.14 MeV 

8.38 




Bs Bg 



0.420 MeV 

3.09 




D* D* 

-^s 



0.156 MeV 

1.15 




Total 



13.6 MeV 

100 




“From PDG Ref. [32]. 

^Using predicted n^Si — n^So splitting and measured mass. 

“3^Pi from LHCb | 33| and 3^P2, 3^Po and 3^Pi using predicted 
splittings with 3®Pi measured mass. 


from the measured values for T{ 2 ^P 2 —>■ 1 ^P 2 + tttt) and 
r(2^Pi —)■ l^Pi + tttt). These partial widths were ob¬ 
tained by first finding the total widths for the Xb 2 i‘ 2 P) 
and xt,i(2P) using the measured BR’s from the PDG [32] 
with our predicted partial widths for El transitions for 
X62 (i)( 2P) ^ 7T(25) and Xb2(i)(2P) ^ 7T(1^). We ob¬ 
tain r(x 62 ) = 122 ± 11 keV and r(xbi) = 62.6 ± 4.0 keV. 
Combining with the measured hadronic BR’s we find 
r(23p2 ^ 13P2 -h tttt) = 0.62 ± 0.12 keV and r(23Pi ^ 
l^Pi -I- tttt) = 0.57 ± 0.09 keV leading to ^2(1,1) = 


1.5 keV and ^o(l) 1) = 1.335 keV. We neglected the small 
differences in phase space between the Xbi and X 62 transi¬ 
tions and remind the reader that model dependence has 
been introduced into these results by using model pre¬ 
dictions for the radiative transition partial widths. For 
example, using the El partial width predictions from 
Kwong and Rosner [S^ results in slightly different to¬ 
tal widths and hadronic transition partial widths. Using 
these values for Aq and A 2 with Eq. [^ we obtain the 
hadronic transition partial widths for 2P —>■ IP transi- 
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TABLE VIII: Partial widths and branching ratios for strong and electromagnetic decays and transitions for the IP states. For 
P-wave annihilation decays M designates R'(0), the first derivative of the radial wavefunction at the origin. See the caption 
to Table [Tv] for further details. 


Initial 

Final 

Mf 

M 

Predicted 

Measured 

state 

state 

(GeV) 


Width (keV) 

BR (%) 

BR (%) 

X62(1®P2) 

99 


1.318 

147 

81.7 


9.912“ 

77 


1.318 

9.3 X 10“® 

5.2 X 10“^ 



T(PS'i)7 

9.460 

1.028 

32.8 

18.2 

19.1 ± 1.2“ 


/ii,(l'Pi)7 

9.899 

1.000 

9.6 X 10“® 

5.3 X 10“® 



Total 



180. 

100 


Xm(1"Pi) 

99 + 9 


1.700 

67 

70. 


9.893“ 

T(CS'i)7 

9.460 

1.040 

29.5 

31 

33.9 ±2.2“ 


Total 



96 

100 


X6o(1^Po) 

99 


2.255 

2.6 MeV 

~ 100 


9.859“ 

77 


2.255 0.15 

5.8 X 10“® 



T(PS'i)7 

9.460 

1.050 

23.8 

0.92 

1.76 ±0.35“ 


Total 



2.6 MeV 

100 


hbil^Pi) 

999 


1.583 

37 

51 


9.899“ 


9.398 

0.922 

35.7 

49 

4g + 8a 


Xh(1^Pi)7 9.893 

1.000 

1.0 X 10“® 

1.4 X 10“® 



Xm(1"Po)7 9.859 

0.998 8.9 X 10“'‘ 

1.2 X 10“® 



Total 



73 

100 



“From PDG Ref. [32). 


tions given in Table IX A note of caution is that Aq and 
A 2 are sensitive to small variations in the input values 
of the partial widths so given the experimental errors on 
the input values, the predictions should only be regarded 
as rough estimates. 

There are no measured BR’s for 3P hadronic tran¬ 
sitions that can be used as input for other 3P transi¬ 
tions. Furthermore, as pointed out by Kuang and Yan 
[60j . hadronic transitions are dependent on intermediate 
states with complicated cancellations contributing to the 
amplitudes so that predictions are rather model depen¬ 
dent. To try to take into account the structure depen¬ 
dence of the amplitudes we make the assumption that 
once phase space and scaling factors (as in Eq. are 
factored out, the ratios of amplitudes will be approxi¬ 
mately the same for transitions between states with the 
same number of nodes in the initial and final states, i.e.: 


nitude estimate of the transitions that can at least tell 
us if the transition is big or small and how significant 
its contribution to the total width will be. With this 
prescription we obtain Ao{3P —>■ IP) = 0.68 keV and 
A2(3P —^ IP) = 4.94 keV. The resulting estimates for 
the 3P —>■ IP hadronic transitions are given in Table 
To obtain these results we used spin averaged P-wave 
masses for the phase space factors given all the other 
uncertainties in these estimates. If we don’t include the 
(r^) rescaling factors, the partial widths increase by 45%, 
which is another reminder that these estimates should be 
regarded as educated guesses that hopefully get the order 
of magnitude right. 


C. IPj —>■ IS -I- TTTT and 2Dj —>■ ID + tttt 


A'jSP ^ IP) A'{3S ^ IS) 

A'(2P ^ IP) "" A'{2S ^ IS) ^ 

where the A'’s have factored out the phase space and 
scaling factors in the amplitude. Thus, we will relate 
the 3P —?► IP partial widths to measured 2P —IP par¬ 
tial widths by rescaling the phase space, the bb separa¬ 
tion factors and using the relationship between ampli¬ 
tudes outlined in Eq. [T^ We understand that this is far 
from rigorous but hope that it captures the gross fea¬ 
tures of the transition and will give us an order of mag¬ 


The BaBar collaboration measured BR{\^D 2 —>■ 

T(IS') -I- 7r+7r“) = (0.66 ± 0.16)% [51]. It should be 
noted that BaBar used the predicted partial widths for 
XbJ' —t from Ref. |3S] as input to obtain this 

value. Combining this measured BR with the remain¬ 
ing 1^P2 partial decay widths given in Table XIV we 


obtain T{1^D2 T)!^) -f Tr+Tr") = 0.169 ± 0.045 keV. 

Eor comparison, using the predictions for the 1^I?2 de¬ 
cays from Ref. [35] we obtain the hadronic width T = 
0.186 ± 0.047 keV. Prior to the BaBar measurement, 
predictions for this transition varied from 0.07 keV to 
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TABLE IX: Partial widths and branching ratios for strong and electromagnetic decays and transitions for the 2P states. For 
P-wave annihilation decays M designates R'(0), the first derivative of the radial wavefunction at the origin. See the caption 
to Table HVl for further details. 


Initial 

Final 

Mf 

M 

Predicted 

Measured 

state 

state 

(GeV) 


Width (keV) 

BR (%) 

BR (%) 

Xb2(2^R2) 

99 


1.528 

207 

89.0 


10.269“ 

77 


1.528 

1.2 X 10"^ 

5.2 X 10"® 



T(23Si)7 

10.023 

1.667 

14.3 

6.15 

10.6 ± 2.6“ 


T(13Si)7 

9.460 

0.224 

8.4 

3.6 

7.0 ±0.7“ 


T3(1"P3)7 

10.172 

-1.684 

1.5 

0.65 



T2(1®P2)7 

10.164 

-1.625 

0.3 

0.1 



Ti(l3Pi)7 

10.154 

-1.561 

0.03 

0.01 



hi,(2iPi)7 

10.260 

1.000 

2.8 X 10"® 

1.2 X 10"® 



hi,(l^Pi)7 

9.899 

-0.011 

2.4 X lO"'^ 

1.0 X lO"'^ 



Xi)2(l^P2)7r7r 



0.62 ±0.12* 

0.27 

0.51 ±0.09“ 


Xi>i(l^-Pi)7r7r 



0.23 

0.10 



Xi>o(l^^o)7i'7r 



0.10 

0.043 



Total 



232.5 

100 



99 + 9 


1.857 

96 

82 


10.255“ 

T(23Si)7 

10.023 

1.749 

13.3 

11.3 

19.9 ± 1.9“ 


T(13Si)7 

9.460 

0.184 

5.5 

4.7 

9.2 ±0.8“ 


T2(1®P2)7 

10.164 

-1.701 

1.2 

1.0 



Ti(l3Pi)7 

10.154 

-1.639 

0.5 

0.4 



hi,(l^Pi)7 

9.899 

-0.035 

2.2 X 10"® 

1.9 X 10"® 



Xi)2(l^P2)7I'7r 



0.38 

0.32 



Xi)i(l^-Pi)7r7r 



0.57 ±0.09* 

0.48 

0.91 ± 0.13“ 


Xi>o(l^Po)7i'7r 



~ 0 

~ 0 



Total 



117 

100 


Xi,o(2^Po) 

99 


2.290 

2.6 MeV 

~ 100 


10.232“ 

77 


2.290 

0.15 

5.7 X 10"® 



T(23Si)7 

10.023 

1.842 

10.9 

0.42 

4.6 ± 2.1“ 


T(13Si)7 

9.460 

0.130 

2.5 

9.6 X 10"^ 

0.9 ±0.6“ 


Ti(l3Pi)7 

10.154 

-1.731 

1.0 

3.8 X 10"^ 



hi,(l^Pi)7 

9.899 

-0.079 

9.7 X 10"® 

3.7 X 10"^ 



Xi)2(l^P2)7r7r 



0.5 

2 X 10"^ 



Xi>i(l^-Pi)^7r 



~ 0 

~ 0 



Xi>o(l^Po)7i'7r 



0.44 

1.7 X 10"^ 



Total 



2.6 MeV 

100 


hb(2^Fi) 

999 


1.758 

54 

64 


10.260“ 

7i,(2^So)7 

9.999 

1.510 

14.1 

17 

48 ± 13“ 


76(l^-S'o)7 

9.398 

0.252 

13.0 

15 

22 ±5“ 


762 ( 1 ^^ 2)7 

10.165 

-1.689 

1.7 

2.0 



X62(1^P2)7 

9.912 

-0.027 

2.2 X 10"® 

2.6 X 10"® 



Xi>i(l^-Pi)7 

9.893 

-0.023 

1.1 X lO"'^ 

1.3 X lO"'^ 



Xi>o(l^Po)7 

9.859 

0.019 

3.2 X lO"'^ 

3.8 X lO"'^ 



/i6(l^Pi)7r7r 



0.94 

1.1 



Total 



84 

100 



“From PDG Ref.|32|. 
'’Input, see text. 
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TABLE X: Partial widths and branching ratios for strong 
and electromagnetic decays and transitions for the 3P states. 
For P-wave annihilation decays M designates P^(0), the first 
derivative of the radial wavefunction at the origin. See the 
caption to Table HYl for further details. 


Initial 

Final 

Mf 

M 

Width 

BR 

state 

state 

(GeV) 


(keV) 

(%) 

Xm(3®P2) 

99 


1.584 

227 

91.9 

10.528'' 

77 


1.584 

1.3 X 10"^ 

5.3 X lO”'^ 


T(33Si)7 

10.355 

2.255 

9.3 

3.8 


T(23Si)7 

10.023 

0.323 

4.5 

1.8 


T(1'^Si)7 

9.460 

0.086 

2.8 

1.1 


T3(23P3)7 

10.455 

-2.568 

1.5 

0.61 


T2(23P2)7 

10.449 

-2.482 

0.32 

0.13 


Ti(23Pi)7 

10.441 

-2.389 

0.027 

0.011 


T3(1"P3)7 

10.172 

0.042 

0.046 

0.019 


Xf)2(l®P2)7r7r 



0.68 

0.28 


Xf>i(l^-Pi)7r7r 



0.52 

0.21 


Xi>o( 1 ^^ 0 ) 71-71 



0.24 

0.10 


Total 



247 

100 

Xi,i(3'^Pi) 

99 + 9 


1.814 

101 

86.3 

10.516“ 

T(3®Si)7 

10.355 

2.388 

8.4 

7.2 


T(23Si)7 

10.023 

0.278 

3.1 

2.6 


T(1'^Si)7 

9.460 

0.061 

1.3 

1.1 


T 2 (23^2)7 

10.449 

-2.595 

1.1 

0.94 


Ti(23Pi)7 

10.441 

-2.506 

0.47 

0.40 


T2(1®P2)7 

10.164 0.060 

0.080 

0.068 


Ti(l3Pi)7 

10.154 0.029 

7.0 X lO"'^ 

6.0 X 10“'^ 


Xf)2(l®P2)7r7r 



0.88 

0.75 


Xi>l(l®-Pl) 71-71 



0.56 

0.48 


Xf>o(l®Po)7171 



~ 0 

~ 0 


Total 



117 

100 

Xf)o(3®Po) 

99 


2.104 

2.2 MeV 

~100 

10.500'' 

77 


2.104 

0.13 

5.9 X 10~'^ 


T(33Si)7 

10.355 

2.537 

6.9 

0.31 


T(23Si)7 

10.023 

0.213 

1.7 

0.077 


T(1'^Si)7 

9.460 

0.029 

0.3 

0.01 


Ti(23Pi)7 

10.441 

-2.637 

1.0 

0.045 


Ti(l3Pi)7 

10.154 0.084 

0.20 

0.0091 


Xi)2(l'^P2)7171 



1.2 

0.054 


Xi>l(l®-Pl) 71-71 



~ 0 

~ 0 


Xf)0(l®Po)7171 



0.27 

0.012 


Total 



2.2 MeV 

100 

hbiS^Pi) 

999 


1.749 

59 

71 

10.519'' 

7?i,(3^So)7 

10.337 

2.047 

8.9 

11 


7i,(2^So)7 

9.999 

0.418 

8.2 

9.9 


76 ( 1 ^ So) 7 

9.398 

0.091 

3.6 

4.3 


762 ( 2 -'^ 2)7 

10.450 

-2.579 

1.6 

1.9 


762(1'P2)7 

10.165 

0.052 

0.081 

0.098 


X62(1®P2)7 

9.912 

-0.032 

1.4 X 10"^ 

1.7 X 10“'' 


X6i(l®Pi)7 

9.893 

-0.031 

9.3 X 10"'' 

1.1 X 10“" 


X6o(l®f’o)7 

9.859 

-0.016 

9.8 X 10"'' 

1.2 X 10“^^ 


/l6(l'Pl)7171 



1.4 

1.7 


Total 



83 

100 


“From LHCb Ref. [33]. 

*’Using the predicted 3P splittings with the measured 3®Pi mass. 


TABLE XL Partial widths and branching ratios for strong 
and electromagnetic decays and transitions and OZI allowed 
strong decays for the 4P states. For P-wave annihilation de¬ 
cays M designates R'{0), the first derivative of the radial 
wavefunction at the origin. Details of the OZI allowed decay 
amplitudes are described in the appendix. See the caption to 
Table HVl for further details. 


Initial 

Final 

Mf 

M 

Width 

BR 

state 

state 

(GeV) 


(keV) 

(%) 

X62(4®P2) 

99 


1.646 

248 

0.569 

10.798 

77 


1.646 

1.5 X 10“" 

3.4 X 10“® 


T(43S'i)7 

10.579 

2.765 

28.1 

6.44 X 10“" 


T(3®5'i)7 

10.355 

0.427 

5.4 

1.2 X 10“" 


T(235'i)7 

10.023 

0.063 

0.59 

1.4 X 10“® 


T(1'^S'i)7 

9.460 

0.056 

2.2 

5.0 X 10“'' 


Ta ( 33 ^ 3)7 

10.711 

-3.310 

4.3 

9.9 X 10“" 


T 2 ( 33 ^ 2)7 

10.705 

-3.202 

0.88 

2.0 X 10“" 


Ti(33Pi)7 

10.698 

-3.084 

0.68 

1.6 X 10“" 


BB 



8.74 MeV 

20.0 


BB* 



28.1 MeV 

64.4 


B*B* 



5.05 MeV 

11.6 


BsBs 



0.593 MeV 

1.36 


BsBt 



0.833 MeV 

1.91 


Total 



43.6 MeV 

100 

X6i(4®A) 

99 + 9 


1.849 

no 

0.36 

10.788 

T(4'^S'i)7 

10.579 

2.942 

27.7 

9.17 X 10“" 


T(3®5'i)7 

10.355 

0.373 

3.8 

1.3 X 10“" 


T(235'i)7 

10.023 

0.035 

0.18 

6.0 X 10“" 


T(1'^S'i)7 

9.460 

0.038 

1.0 

3.3 X 10“" 


T 2 ( 33 ^ 2)7 

10.705 

-3.345 

3.4 

1.1 X 10“" 


Ti(33Pi)7 

10.698 

-3.234 

1.4 

4.6 X 10“" 


BB* 



20.6 MeV 

68.3 


B*B* 



0.478 MeV 

1.58 


BsBl 



8.93 MeV 

29.6 


Total 



30.2 MeV 

100 

X6o(4®Po) 

99 


2.079 

2.1 MeV 

6.1 

10.775 

77 


2.079 

0.13 

3.8 X 10“" 


T(43S'i)7 

10.579 

3.139 

26.0 

7.54 X 10“" 


T(33S'i)7 

10.355 

0.295 

2.2 

6.4 X 10“" 


T(23S'i)7 

10.023 

-0.001 

9 X 10“® 

3 X 10“'" 


T(13S'i)7 

9.460 

0.017 

0.21 

6.1 X 10“" 


Ti(33Pi)7 

10.698 

-3.397 3.8 

1.1 X 10“" 


BB 



20.0 MeV 

58.0 


B*B* 



12.2 MeV 

35.4 


BsBs 



0.129 MeV 0.374 


Total 



34.5 MeV 

100 

h6(4'Pi) 

999 


1.790 

64 

0.16 

10.790 

76(4'So)7 

10.567 2.808 

24.4 

6.07 X 10“" 


?76(3'5'o)7 

10.337 0.587 

10.8 

2.69 X 10“" 


76(2'So)7 

9.999 

0.055 

0.48 

1.2 X 10“" 


76(l'So)7 

9.398 

0.052 

2.1 

5.2 X 10“" 


762(3'P2)7 10.790 

-3.325 4.7 

1.2 X 10“" 


BB* 



31.8 MeV 

79.1 


B*B* 



4.09 MeV 

10.2 


BsBl 



4.18 MeV 

10.4 


Total 



40.2 MeV 

100 
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TABLE XII: Partial widths and branching ratios for strong 
and electromagnetic transitions and OZI allowed strong de¬ 
cays for the 5®P2 and 5®Pi states. Details of the OZI allowed 
decay amplitudes are described in the appendix. See the cap¬ 
tion to Table HVl for further details. 


TABLE XIII: Partial widths and branching ratios for strong 
and electromagnetic transitions and OZI allowed strong de¬ 
cays for the 5®Po and 5^Pi states. Details of the OZI allowed 
decay amplitudes are described in the appendix. See the cap¬ 
tion to Table nYi for further details. 


Initial 

state 

Final 

state 

Mf 

(GeV) 

M 

Width 

(keV) 

BR 

(%) 

Initial 

state 

Final 

state 

Mf 

(GeV) 

M 

Width 

(keV) 

BR 

(%) 

Xf>2(5®P2) 

T(5®Si)7 

10.876 

3.232 

11.5 

2.06 X 10-2 

Xi)o(52Po) 

T(52Si)7 

10.876 

3.668 

10.0 

1.86 X 10-2 

11.022 

T(43Si)7 

10.579 

0.595 

10.4 

1.86 X 10-2 

11.004 

T(42S'i)7 

10.579 

0.446 

5.2 

9.6 X 10-2 


T(3®Si)7 

10.355 

0.024 

0.06 

1 X 10-"^ 


T(32Si)7 

10.355 

-0.043 

0.16 

3.0 X 10-'‘ 


T(23Si)7 

10.023 

0.067 

1.4 

2.5 X 10-2 


T(22S'i)7 

10.023 

0.035 

0.36 

6.7 X 10-'‘ 


T(USi)7 

9.460 

0.042 

1.9 

3.4 X 10-2 


T(12Si)7 

9.460 

0.014 

0.22 

4.1 X lO-'^ 


T3(43D3)7 

10.939 

-3.955 

5.4 

9.7 X 10-2 


Ti(42Di)7 

10.928 

-4.038 

5.1 

9.5 X 10-2 


T2(43D2)7 

10.934 

-3.828 

1.1 

2.0 X 10-2 


BB 



4.52 MeV 

8.38 


Ti(43Di)7 

10.928 

-3.685 

0.08 

1 X 10-"^ 


BP(lPi) 



7.16 MeV 

13.3 


BB 



0.456 MeV 

0.816 


B*B* 



40.3 MeV 

74.8 


BB* 



2.71 MeV 

4.85 


Bs Bs 



0.166 MeV 0.308 


BB(lPi) 



4.72 MeV 

8.44 


BtBt 



1.71 MeV 

3.17 


BB{1P{) 



15.8 MeV 

28.3 


Total 



53.9 MeV 

100 


B*B* 



31.3 MeV 

56.0 

hb{5^Pi) 

rjbi^^Soh 

10.867 2.900 

9.8 

2.1 X 10-2 


BsBs 



0.154 MeV 

0.275 

11.016 

?7i)(4^So)7 

10.567 0.824 

20.8 

4.54 X 10-2 


BsBt 



0.130 MeV 

0.232 


Vbi^^Soh 

10.337 

-0.003 

9.7 X 10-^ 

2.1 X 10-2 


BlBt 



0.618 MeV 

1.10 


Vbi2^Soh 

9.999 

0.030 

0.29 

6.3 X 10-'‘ 


Total 



55.9 MeV 

100 


?7i)(l^So)7 

9.398 

0.042 

2.2 

4.8 X 10-2 

Xm(5®Pi) 

T(53Si)7 

10.876 

3.439 

11.0 

1.75 X 10-2 


77i,2 ( 42 ^ 2)7 10.790 

-3.967 6.0 

1.3 X 10-2 

11.014 

T(43Si)7 

10.579 

0.534 

8.0 

1.3 X 10-2 


BB* 



11.7 MeV 

25.5 


T(3®Si)7 

10.355 

-0.006 

2.8 X 10"^ 

4.4 X 10-® 


BB(l2Po) 



6.81 MeV 

14.9 


T(23Si)7 

10.023 

0.052 

0.83 

1.3 X 10-2 


BB(lPi) 



0.412 

9.00 X 10-'‘ 


T(USi)7 

9.460 

0.029 

0.90 

1.4 X 10-2 


BB{1P[) 



0.689 

1.50 X 10-2 


T 2 ( 43 ^ 2)7 

10.934 

-3.990 

4.4 

7.0 X 10-2 


B*B* 



26.4 MeV 

57.6 


Ti(43Di)7 

10.928 

-3.857 

1.7 

2.7 X 10-2 


BsBl 



55.8 

0.122 


BB* 



16.7 MeV 

26.5 


BlBl 



0.670 MeV 

1.46 


BB{^Pg) 



0.306 

4.86 X 10-^ 


Total 



45.8 MeV 

100 


BB(lPi) 



13.5 MeV 

21.4 








BB{1P[) 



6.82 MeV 

10.8 








B*B* 



25.1 MeV 

39.8 








BsBt 

BlBl 

Total 


3.41 X 10~ 


21.5 

0.614 MeV 0.975 
63.0 MeV 100 


24 keV [Sni [Ml EHl m] • Using the partial width value F = 
0.169±0.045 keV as input we obtain ^ 2 ( 2 , 0) = 0.845 keV 
which we use along with phase space (Eq. and (r^) 
rescaling factors (Eqj8|to obtain the ID —>■ IS'tttt transi¬ 
tions given in Table jXIV For comparison we include the 
measurements for the l^Di and l^Ha transitions from 
Ref. [M] which are less certain than those for the 1^D2 
transitions. 

Because we have no data on the 2D states we use the 
same strategy to estimate 2D ID transitions as we 


did in estimating the 3P —>■ IP transitions. We assume 
that rescaling an amplitude with the same number of 
nodes in the initial and final state wavefunctions will cap¬ 
ture the gross features of the complicated overlap inte¬ 
grals with intermediate wavefunctions. We use the Aq 
and A 2 amplitudes from the 2P —IP transitions as 
input and rescale the amplitudes using the appropriate 
phase space factors and (r^) rescaling factors. This gives 
Ao(2, 2) = 3.1 X 10-2 keV and ^ 2 ( 2 ,2) = 4.6 x lO"® keV. 
As a check we also estimated Ag found by rescaling the Ag 
amplitude obtained from the 2S —> IS' transition where 
only Ag contributes and found it to be roughly a factor 
of 40 smaller than the value obtained from the 2P —IP 
transition. This should be kept in mind when assessing 
the reliability of our predictions. In any case, the partial 
widths obtained for the 2D —)■ ID hadronic transitions 


are sufficiently small (see Table XV) that the large un- 
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TABLE XIV: Partial widths and branching ratios for strong and electromagnetic decays and transitions for the ID states. For 
D-wa.ve annihilation decays M designates R"{0), the second derivative of the radial wavefunction at the origin. See the caption 
to Table [Tv] for further details. 


Initial 

Final 

Mf 

M 

Predicted 

Measured 

state 

state 

(GeV) 


Width (keV) 

BR (%) 

BR (%) 

T3(1"F3) 

Xi>2(l®F2)7 

9.912 

1.830 

24.3 

91.0 


10.172*’ 

?7i,2(l^F2)7 

10.165 

1.000 

1.5 X 10"® 

5.6 X 10“® 



ggg 


0.9923 

2.07 

7.75 



T(F5'i)7r+7r- 



0.197 

0.738 

0.29 ± 0.23 (or < 0.62 at 90% C.L.)“ 


Total 



26.7 

100 


T2(1*'F2) 

Xi>2(l®F2)7 

9.912 

1.835 

5.6 

22 


10.164“ 

Xi>i(l®Fi)7 

9.893 

1.762 

19.2 

74.7 



ggg 


1.149 

0.69 

2.7 



T(G5'i)7r+7r- 



0.169 ±0.045 

0.658 

0.66 ±0.16“ 


Total 



25.7 

100 


Ti(FFi) 

i+r 



1.38 eV 

3.93 X 10"® 


lO.lSS*" 

Xi>2(l®F2)7 

9.912 

1.839 

0.56 

1.6 



Xi>i(l®Fi)7 

9.893 

1.768 

9.7 

28 



Xi>o(l®Fo)7 

9.859 

1.673 

16.5 

47.1 



ggg 

1.356 


8.11 

23.1 



T(l35'i)7r+7r- 



0.140 

0.399 

0.42 ± 0.29 (or < 0.82 at 90% C.L.)“ 


Total 



35.1 

100 


77i,2(l^F2) 

/ii,(l^Fi)7 

9.899 

0.178 

24.9 

91.5 


lo.ies*" 

r7i,(l*-So)7r+7r“ 



0.35 

1.3 



gg 


1.130 

1.8 

6.6 



Total 



27.2 

100 



“From BaBar |84| . 

*’Using predicted splittings and DD 2 mass from Ref. 1841 . 
“See Section IV C for the details of how this was obtained. 


certainties will not change our conclusions regarding the 
2D states. 


D. IF —>■ IP + TTTT and IG —t ID + tttt 

We take the same approach as we used to estimate 
some of the hadronic transition widths given above. 
We take a measured width, in this case the 1^D2 —t 
T(lS')7r+7r“, and rescale it using ratios of phase space 
factors and separation factors to estimate the IF —>■ IF 
and 1C? —>■ ID transitions. We obtain A2{1F —>■ IF) = 
0.027 keV and A 2 {IG ID) = 1.94 x IQ-^ keV. These 
small values are primarily due to the ratio of phase space 
factors which roughly go like the mass difference to the 
7th power. Given that the IF—IF, IF—IF and 1C? —IF 
mass splittings are ^ 0.69, 0.46 and 0.38 MeV respec¬ 
tively resulting in little available phase space one can un¬ 
derstand why the amplitudes are small. We include our 
estimates for these transitions in Tables IXIXI and I XXIII 
While the estimates may be crude the point is that we 
expect these partial widths to be quite small. 


VI. STRONG DECAYS 

For states above the BB threshold, we calculate OZI 
allowed strong decay widths using the ^Fq quark pair 
creation model [laiisiiiiiiiTiiHS] which proceeds through 
the production of a light qq pair {q = u, d, s) followed 
by separation into BB mesons. The qq pair is assumed 
to be produced with vacuum quantum numbers (0++). 
There are a number of predictions for T strong decay 
widths in the literature using the ^Fq model [HI HZ] and 
other models [5S] , but a complete analysis of their strong 
decays had yet to be carried out prior to this work. We 
give details regarding the notation and conventions used 
in our ^Fq model calculations in Appendix [A| to make it 
more transparent for an interested reader to reproduce 
our results. 

We use the meson masses listed in Tables Ul lIIIl If avail¬ 
able, the measured value, M^xp, is used as input for cal¬ 
culating the strong decay widths, rather than the pre¬ 
dicted value, Mtheo- When the mass of only one meson 
in a multiplet has been measured, we estimate the input 
masses for the remaining states following the procedure 
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TABLE XV: Partial widths and branching ratios for strong 
and electromagnetic decays and transitions for the 2D states. 
For D-wave annihilation decays M designates R"(0), the sec¬ 
ond derivative of the radial wavefunction at the origin. See 
the caption to Table m for further details. 


Initial 

Final 

Mf 

M 

Width 

BR 

state 

state 

(GeV) 


(keV) 

(%) 

T3(23D3) 

Xf>2(2®P2)7 

10.269 

2.445 

16.4 

65.1 

10.455 

Xm(1®P2)7 

9.912 

0.200 

2.6 

10. 


Xf>4(l®P4)7 

10.358 

-1.798 

1.7 

6.7 


Xi)3(l^P3)7 

10.355 

-1.751 

0.16 

0.63 


Xi>2(l®P2)7 

10.350 

-1.702 

5 X 10“® 

2 X 10"^ 


762(2^1)2)7 

10.450 0.999 

6.5 X 10"® 

2.6 X 10“® 


762(1^1)2)7 

10.165 

-0.033 

1.1 X 10“^ 

4.4 X 10“® 


ggg 


1.389 

4.3 

17 


T3{l^D3)mv 



7.4 X 10"® 

2.9 X 10“^ 


T2(l®D2)7r7r 



1.9 X 10"® 

7.5 X 10“® 


Ti(l®Di)7r7r 



1.9 X 10"'^ 

7.5 X 10“^ 


Total 



25.2 

100 

T2{2^D2) 

X62(2®P2)7 

10.269 

2.490 

3.8 

17 

10.449 

Xbi(2®Pi)7 

10.255 

2.359 

12.7 

56.2 


X62(1®P2)7 

9.912 

0.161 

0.4 

2 


Xbi(l®-Pi)7 

9.893 

0.224 

2.6 

12 


X63(1^P3)7 

10.355 

-1.806 

1.5 

6.6 


Xb2(l®P2)7 

10.350 

-1.758 0.21 

0.93 


762(1^^2)7 

10.165 

-0.047 2.1 X 10"® 

9.3 X 10“® 


ggg 


1.568 

1.4 

6.2 


T3(l®P3)7r7r 



2.6 X 10"® 

1.2 X 10“® 


T2(l®P2)7r7r 



7.4 X 10"® 

3.3 X 10“2 


Ti(l®Pi)7r7r 



2.3 X 10"® 

1.0 X 10“® 


Total 



22.6 

100 

Ti(23Di) 

i+r 



1.99 eV 

5.28 X 10"® 

10.441 

X62(2®P2)7 

10.269 

2.535 

0.4 

1 


X6i(2®Pi)7 

10.255 

2.409 

6.5 

17 


X6o(2®Po)7 

10.232 

2.243 

10.6 

28.1 


X62(1®P2)7 

9.912 

0.118 

0.02 

0.05 


X6i(l®Pi)7 

9.893 

0.184 

0.9 

2 


X6o(l®Po)7 

9.859 

0.260 

2.9 

7.7 


X62(1®P2)7 

10.350 

-1.815 

1.6 

4.2 


762(1^^2)7 

10.165 

-0.061 

CO 

CO 

X 

0 

1 

W 

8.8 X 10“® 


ggg 


1.771 

14.8 

39.2 


T3(l®P3)7r7r 



4.4 X 10"'^ 

1.2 X 10“® 


T2(l®P2)7r7r 



3.9 X 10“® 

1.0 X 10“® 


Ti(l3p>i)7r7r 



7.4 X 10"® 

2.0 X 10“2 


Total 



37.7 

100 

r]h2{2^ D 2 ) 

h6(2iPi)7 

10.260 

2.390 

16.5 

67.1 

10.450 

h6(l^Pi)7 

9.899 

0.212 

3.0 

12 


ft63(l^P3)7 

10.355 

-1.802 

1.8 

7.3 


T3(1"P3)7 

10.172 

-0.072 

6.5 X 10"® 

2.6 X 10“^ 


T2(1®P2)7 

10.164 

-0.048 

2.2 X 10“® 

8.9 X 10“® 


Ti(l3Pi)7 

10.154 

-0.046 

1.4 X 10"® 

5.7 X 10“® 


gg 


1.530 

3.3 

13 


762(l^P2)7r7r 



7.4 X 10"® 

3.0 X 10“® 


Total 



24.6 

100 


TABLE XVI: Partial widths and branching ratios for strong 
and electromagnetic decays and transitions and OZI allowed 
decays for the 3D states. For D-wave annihilation decays M 
designates R"{0), the second derivative of the radial wave- 
function at the origin. Details of the OZI allowed decays 
amplitudes are described in the appendix. See the caption to 
Table HVl for further details. 


Initial 

Final Mf A4 

Width 

BR 

state 

state (GeV) 

(keV) 

(%) 

T3(3®P3) 

Xf>2(3®P2)7 10.528 3.022 

23.6 

1.19 X 10"® 

10.711 

Xi>2(2®P2)7 10.269 0.265 

2.5 

1.3 X 10"® 


Xi>2(l®P2)7 9-912 0.064 

0.82 

4.1 X lO"'^ 


Xm(2®P4)7 10.622 -2.683 3.0 

1.5 X 10"® 


Xi>3(2®P3)7 10.619 -2.617 0.27 

1.4 X lO""' 


Xi>2(2®P2)7 10.615 -2.545 

8 X 10“® 

4 X 10“® 


ggg i-69i 

6.6 

3.3 X 10"® 


BB 

16.3 MeV 

8.23 


BB* 

72.9 MeV 

36.8 


B*B* 

109 MeV 

55.0 


Total 

198 MeV 

100 

T2(3®P2) 

Xf>2(3®P2)7 10.528 3.098 

5.6 

4.3 X 10"® 

10.705 

Xm(3®Pi)7 10.516 2.919 

18.2 

1.41 X 10"® 


Xi>2(2®P2)7 10.269 0.218 

0.40 

3.1 X lO"'^ 


Xm(2®Pi)7 10.255 0.303 

2.5 

1.9 X 10"® 


Xi>2(l®P2)7 9.912 0.043 

0.09 

7 X 10~® 


Xbiil^Pih 9.893 0.062 

0.59 

4.6 X lO"'^ 


Xi>3(2®P3)7 10.619 -2.698 

2.6 

2.0 X 10"® 


Xi>2(2®P2)7 10.615 -2.628 0.36 

2.8 X lO""' 


ggg 1-875 

2.0 

1.6 X 10"® 


BB* 

52.4 MeV 

40.6 


B*B* 

76.5 MeV 

59.3 


Total 

128.9 MeV 

100 

Ti(3®Pi) 

i+r 

2.38 eV 

2.30 X 10"® 

10.698 

Xi>2(3®P2)7 10-528 3.174 

0.58 

5.6 X lO"'^ 


Xm(3®Pi)7 10.516 3.003 

9.5 

9.2 X 10"® 


Xi>o(3®Po)7 10.500 2.775 

14.0 

1.35 X 10"® 


Xi>2(2®P2)7 10.269 0.165 

0.02 

2 X 10“® 


Xm(2®Pi)7 10.255 0.256 

0.96 

9.3 X 10"^ 


Xw(2®Po)7 10.233 0.354 

2.8 

2.7 X 10"® 


Xw(l®Pi)7 9.893 0.040 

0.13 

1.3 X lO"'^ 


Xf>o(l®Po)7 9.859 0.069 

0.59 

5.7 X lO"'^ 


Xi>2(2®P2)7 10.615 -2.712 

2.7 

2.6 X 10"® 


Xf>2(l®P2)7 10.350 0.039 

0.39 

3.8 X 10"® 


ggg 2.O8I 

21.2 

2.05 X 10"® 


BB 

23.8 MeV 

23.0 


BB* 

0.245 MeV 0.236 


B*B* 

79.5 MeV 

76.7 


Total 

103.6 MeV 

100 

rib2{3^D2) 

hi,i(3®Pi)7 10.519 2.956 

24.1 

1.43 X 10"® 

10.706 

/ii,i(2®Pi)7 10.260 0.285 

2.9 

1.7 X 10"® 


/ii,i(l®Pi)7 9.899 0.061 

0.76 

4.5 X lO"'^ 


hi,3(2®P3)7 10.619 -2.691 

3.1 

1.8 X 10"® 


/ii,2(l®P2)7 10.355 0.030 

0.02 

1 X 10~® 


gg 1-839 

4.7 

2.8 X 10"® 


BB* 

77.8 MeV 

46.1 


B*B* 

90.9 MeV 

53.9 


Total 

168.7 MeV 

100 
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TABLE XVII: Predicted partial widths and branching ratios 
for strong and electromagnetic transitions and OZI allowed 
decays for the 4^Da and 4^_D2 states. Details of the OZI 
allowed decay amplitudes are described in the appendix. See 
the caption to Table HYl for further details. 


TABLE XVIII: Predicted partial widths and branching ra¬ 
tios for strong and electromagnetic transitions for 4^Di and 
4^D2 states. Details of the OZI allowed decay amplitudes are 
described in the appendix. See the caption to Table |IV| for 
further details. 


Initial 

state 

Final 

state 

Mf 

(GeV) 

M 

Width 

(keV) 

BR 

(%) 

Initial 

state 

Final 

state 

Mf 

(GeV) 

M 

Width 

(keV) 

BR 

(%) 

T3(43D3) 

Xf>2(4®P2)7 

10.798 

3.553 

15.0 

2.59 X 10“^ 

Ti(43Di) 

i+r 



2.18 eV 

3.04 X 10"® 

10.939 

Xi>2(3®P2)7 

10.528 

0.322 

2.9 

5.0 X 10“® 

10.928 

Xi)2(4®P2)7 

10.798 

3.756 

0.36 

5.0 X lO"'^ 


Xf>2(2®P2)7 

10.269 

0.073 

0.63 

1.1 X 10"® 


Xm(4®Pi)7 

10.788 

3.538 

6.1 

8.5 X 10"® 


Xf>2(l®P2)7 

9.912 

0.035 

0.50 

8.6 X 10"^ 


Xw(4®Po)7 

10.775 

3.262 

9.0 

1.2 X 10"^ 


Xi,4(3®F4)7 

10.856 

-3.143 

3.9 

6.7 X 10“^ 


Xi>2(3®P2)7 

10.528 

0.210 

0.03 

4 X 10~® 


Xf)3 ( 3 ^ 73)7 

10.853 

-3.330 

0.36 

6.2 X 10"^ 


Xm(3®Pi)7 

10.516 

0.331 

1.3 

1.8 X 10"® 


Xf>2(3®7’2)7 

10.850 

-3.241 

0.01 

2 X 10“® 


Xw(3®Po)7 

10.500 

0.476 

4.0 

5.6 X 10"® 


BB 



0.726 MeV 

1.25 


Xi>2(2®P2)7 

10.269 

0.022 

1.5 X 10-® 

2.1 X 10"® 


BB* 



2.94 MeV 

5.07 


Xm(2®Pi)7 

10.255 

0.037 

0.07 

1 X 10“" 


B*B* 



51.5 MeV 

88.8 


Xi,o(2®Po)7 

10.233 

0.060 

0.26 

3.6 X 10“" 


BsBs 



0.265 MeV 

0.457 


Xm(1®Ri)7 

9.893 

0.033 

0.19 

2.6 X 10“" 


BsBl 



0.0827 MeV 0.142 


Xw(l®Po)7 

9.859 

0.054 

0.75 

l.O X 10“® 


BlBl 



2.44 MeV 

4.21 


Xi)2(3®P2)7 

10.850 

-3.448 

3.6 

5.0 X 10“® 


Total 



58.0 MeV 

100 


BB 



3.85 MeV 

5.36 

T2(43D2) 

Xf>2(4®P2)7 

10.798 

3.655 

3.6 

5.6 X 10“® 


BB* 



14.0 MeV 

19.5 

10.934 

Xi>i(4®Ri)7 

10.788 

3.428 

11.6 

1.80 X 10“^ 


B*B* 



50.6 MeV 

70.5 


Xf>2(3®P2)7 

10.528 

0.270 

0.50 

7.8 X 10"^ 


BsBs 



O.lOl MeV 0.141 


Xi>i(3®Pi)7 

10.516 

0.384 

3.3 

5.1 X 10“^ 


BsBt 



0.332 MeV 0.462 


Xi>2(2®P2)7 

10.269 

0.048 

0.07 

1 X 10-^ 


BtBt 



2.94 MeV 

4.09 


Xi,i(2®Pi)7 

10.255 

0.061 

0.35 

5.4 X 10"^ 


Total 



71.8 MeV 

100 


Xi>2(l®P2)7 

9.912 

0.022 

0.05 

8 X 10“® 

rjb2{3} D 2 ) 

hb{4}P^_)'y 

10.790 

3.477 

15.6 

2.58 X 10“^ 


Xi>i(l®Ri)7 

9.893 

0.046 

0.68 

1.1 X 10"® 

10.935 


10.519 

0.362 

3.8 

6.3 X 10“® 



10.853 

-3.431 

3.6 

5.6 X 10"® 


hb{2^Pih 

10.260 

0.065 

0.51 

8.4 X 10“" 


Xi)2(3®7’2)7 

10.850 

-3.344 

0.47 

7.3 X 10“^ 


hbil^Pih 

9.899 

0.042 

0.75 

1.2 X 10“® 


BB* 



25.7 MeV 

40.0 


hi,3(3^P3)7 

10.853 

-3.423 

4.1 

6.8 X 10“® 


B*B* 



36.4 MeV 

56.6 


BB* 



19.4 MeV 

32.1 


BsBl 



0.357 MeV 

0.56 


B*B* 



38.9 MeV 

64.3 


BlBl 



1.80 MeV 

2.80 


BsBl 



0.239 MeV 

0.395 


Total 



64.3 MeV 

100 


BlBl 



1.92 MeV 

3.17 








Total 



60.5 MeV 

100 


described at the end of Sec. m 

We use simple harmonic oscillator wave functions with 
the effective oscillator parameter, j3, obtained by equat¬ 
ing the rms radius of the harmonic oscillator wavefunc- 
tion for the specified (n,/) quantum numbers to the rms 
radius of the wavefunctions calculated using the rela¬ 
tivized quark model of Ref. m- The effective harmonic 
oscillator wavefunction parameters found in this way are 
listed in the final column of Tables mimi For the con¬ 
stituent quark masses in our calculations of both the 
meson masses and of the strong decay widths, we use 
rub = 4.977 GeV, m, = 0.419 GeV, and to, = 0.220 GeV 
{q = u,d). Finally, we use “relativistic phase space” as 


described in Ref. [HI [SS] and in Appendix 

Typical values of the parameters /3 and 7 are found 
from fits to light meson decays [MiMiinn]. The predicted 
widths are fairly insensitive to the precise values used for 
(3 provided 7 is appropriately rescaled. However 7 can 
vary as much as 30% and still give reasonable overall fits 
of light meson decay widths [ESI [SO]. This can result in 
factor of two changes to predicted widths, both smaller 
or larger. In our calculations of Dg meson strong decay 
widths in [18], we used a value of 7 = 0.4, which has also 
been found to give a good description of strong decays 
of charmonium m- However, we found that this value 
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TABLE XIX: Predicted partial widths and branching ratios 
for strong and electromagnetic decays and transitions for IE- 
wave states. For E-wave annihilation decays M designates 
R'"(0), the third derivative of the radial wavefunction at the 
origin. See the caption to Table m for further details. 


Initial 

Final 

Mf 

M 

Width 

BR 

state 

state 

(GeV) 


(keV) 

(%) 

Xi)4(l^E4) 

T3(1"E3)7 

10.172 

2.479 

18.0 

~ 100 

10.358 

XM(l®E2)7r7r 



3.9 X 10"® 

2.2 X 10“^ 


99 


0.868 0.048 

0.27 


Total 



18.0 

100 

Xm(1'F3) 

T3(1"E3)7 

10.172 

2.482 

1.9 

10. 

10.355 

T2(13E2)7 

10.164 

2.442 

16.7 

89.3 


99 


0.974 0.060 

0.32 


Xf)2(l®E2)7r7r 



1.3 X 10"® 

7.0 X 10“^ 


Xf>i(l®7’i)7r7r 



2.6 X 10"® 

1.4 X 10“^ 


Total 



18.7 

100 

Xi>2(l®E2) 

T3(1"E3)7 

10.172 

2.485 

0.070 

0.35 

10.350 

T2(1®E2)7 

10.164 

2.446 

2.7 

14 


Ti(13Ei)7 

10.154 

2.402 

16.4 

82.4 


99 


1.091 

0.70 

3.5 


Xi)2(l®E2)7r7r 



2.6 X 10"'‘ 

1.3 X 10“® 


Xi>i(l®-Pi)7r7r 



1.8 X 10"® 

9.0 X 10“^ 


Xf)o(l®Po)7i'7r 



1.8 X 10"® 

9.0 X 10“^ 


Total 



19.9 

100 

h,3(l^F3) 

7i)2(l^E2)7 

10.165 

2.449 

18.8 

~100 

10.355 

/ii,(l^Ei)7r7r 



3.9 X 10"® 

2.1 X 10“^ 


Total 



18.8 

100 


underestimated the bottomonium strong decay widths 
when compared to the PDG values for the T(4S'), T(5S') 
and T(65') widths. Therefore, we used a value of 7 = 
0.6 in our strong decay width calculations in this paper, 
which was determined by fitting our results to the PDG 
values in the T sector. This scaling of the value of 7 
in different quarkonia sectors has been studied in m- 


The resulting strong decay widths are listed in Tables VI 


IXXIIIl in which we use a more concise notation where BB 
refers to the BB decay mode, BB* refers to BB* +BB*, 
etc. 


We note that our results differ from the recent work 
of Ferretti and Santopinto [5S] , in some cases quite sub¬ 
stantially. This is primarily due to the values chosen 
for the harmonic oscillator parameter fi (with a corre¬ 
sponding change in the pair creation strength 7 ). In 
our calculations we used a value for /? found by fitting 
the rms radius of a harmonic oscillator wavefunction to 
the “exact” wavefunction for each state while Ferretti 
and Santopinto used a common value for all states. An¬ 
other reason our results differ is because Ferretti and San¬ 
topinto included an additional Gaussian smearing func¬ 


TABLE XX: Predicted partial widths and branching ratios for 
strong decays and electromagnetic transitions for the 2E-wave 
states. For E-wave annihilation decays M designates 
the third derivative of the radial wavefunction at the origin. 
Details of the OZI allowed decay amplitudes are described in 
the appendix. See the caption to Table [IV| for further details. 


Initial 

Final 

Mf 

M 

Width 

BR 

state 

state 

(GeV) 


(keV) 

(%) 

Xm(23E4) 

T3(23E3)7 

10.455 

3.053 

19.6 

0.700 

10.622 

T3(1"E>3)7 

10.172 

0.191 

1.4 

5.0 X 10"2 


T5(1"G5)7 

10.532 

-1.886 

1.5 

5.4 X 10"^ 


T4(1"G4)7 

10.531 

-1.848 

0.08 

3 X 10"® 


T3(1"G3)7 

10.529 

-1.808 

1 X 10"® 

4 X 10"® 


99 


1.455 

0.13 

4.6 X 10"® 


BB 



2.73 MeV 

97.5 


BB* 



0.0462 MeV 

1.65 


Total 



2.80 MeV 

100 

Xi>3(2®E3) 

T3(23e>3)7 

10.455 

3.084 

2.1 

1.4 X 10"^ 

10.619 

T2(23E)2)7 

10.449 

3.009 

17.9 

0.116 


T3(1®E3)7 

10.172 

0.156 

0.1 

6 X 10"'^ 


T2(1®E>2)7 

10.164 

0.199 

1.4 

9.1 X 10"® 


T4(1®G4)7 

10.531 

-1.892 

1.4 

9.1 X 10"® 


T3(1"G3)7 

10.529 

-1.852 

0.10 

6.5 X 10-* 


99 


1.583 

0.16 

1.0 X 10"® 


BB* 



15.4 MeV 

~ 100 


Total 



15.4 MeV 

100 

Xi>2(2®E2) 

r3{2^D3h 

10.455 

3.114 

0.08 

9 X 10"® 

10.615 

T2(23E>2)7 

10.449 

3.042 

3.0 

3.4 X 10"® 


Ti(23e>i)7 

10.441 

2.961 

17.5 

1.98 X 10"® 


T3(1"E»3)7 

10.172 

0.118 

2 X 10"® 

2 X 10"® 


T2(1"E>2)7 

10.164 

0.163 

0.16 

1.8 X lO"'^ 


Ti(13e>i)7 

10.154 

0.210 

1.6 

1.8 X 10"® 


T3(1"G3)7 

10.529 

-1.898 

1.4 

1.6 X 10"® 


99 


1.752 

1.77 

2.00 X 10"® 


BB 



83.4 MeV 

94.1 


BB* 



5.20 MeV 

5.987 


Total 



88.6 MeV 

100 

hi,3(2^E3) 

762(2^^2)7 

10.450 

3.019 

19.9 

0.169 

10.619 

762 ( 1 ^ 772)7 

10.148 

0.196 

1.6 

1.4 X 10"® 


7i,4(l^G4)7 

10.530 

-1.890 

1.5 

1.3 X 10"® 


BB* 



11.8 MeV 

~ 100 


Total 



11.8 MeV 

100 


tion in their momentum-space wavefunction overlap to 
model the non-point-like nature of the created qq pair. As 
a numerical check of our programs we reproduced their 
results using their parameters and including the Gaus¬ 
sian smearing function, although we found that the latter 
had little effect on our results. We believe our approach 
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TABLE XXL Partial widths and branching ratios for electro¬ 
magnetic transitions and OZI allowed decays for the 3-F-wave 
states. Details of the OZI allowed decay amplitudes are de¬ 
scribed in the appendix. See the caption to Table |IV| for 
further details. 


Initial 

Final 

Mf 

M 

Width 

BR 

state 

state 

(GeV) 


(keV) 

(%) 

Xbi{S^ F 4 ) 

T3(23D3)7 

10.711 

3.593 

17.9 

1.87 X 10"^ 

10.856 

T3(23D3)7 

10.455 

0.256 

1.9 

2.0 X 10“® 


T3(1"D3)7 

10.172 

0.050 

0.34 

3.6 X 10-3 


T5(23G5)7 

10.772 

-2.776 

2.6 

2.7 X 10“® 


T4(23G4)7 

10.771 

-2.722 

0.14 

1.5 X 10“3 


T3(23G3)7 

10.769 

-2.664 

2.2 X 10“3 

2.3 X 10“® 


BB 



2.84 MeV 

2.97 


BB* 



0.681 MeV 

0.713 


B*B* 



85.7 MeV 

89.7 


BsBs 



0.733 MeV 

0.768 


BsBl 



1.14 MeV 

1.19 


BlBt 



4.43 MeV 

4.64 


Total 



95.5 MeV 

100 

XbsiS^Fa) 

T3(23D3)7 

10.711 

3.646 

1.9 

1.9 X 10“® 

10.853 

T2(23D2)7 

10.705 

3.542 

16.4 

1.62 X 10"^ 


T3(23D3)7 

10.455 

0.214 

0.14 

1.4 X 10“3 


T 2 ( 23772)7 

10.449 

0.271 

1.9 

1.9 X 10“® 


T4(23G4)7 

10.771 

-2.785 

2.4 

2.4 X 10“® 


T3(23G3)7 

10.709 

-2.728 

0.17 

1.7 X 10“3 


BB* 



43.8 MeV 

43.2 


B*B* 



52.4 MeV 

51.7 


BsBt 



3.83 MeV 

3.78 


BlBt 



1.30 MeV 

1.28 


Total 



101.4 MeV 

100 

Xb2i3^F2) 

T3(33773)7 

10.711 

3.699 

0.07 

6 X 10"® 

10.850 

T2(33772)7 

10.705 

3.598 

2.8 

2.6 X 10“® 


Ti(3377i)7 

10.698 

3..487 

16.3 

1.50 X 10"^ 


Ti(2377i)7 

10.441 

0.290 

2.1 

1.9 X 10“® 


T3(23G3)7 

10.769 

-2.794 2.5 

2.3 X 10“® 


BB 



7.85 MeV 

7.21 


BB* 



32.0 MeV 

29.4 


B*B* 



66.0 MeV 

60.6 


BsBs 



0.709 

6.51 X 10"3 


BsBl 



2.50 MeV 

2.30 


BlBt 



0.557 MeV 0.511 


Total 



108.9 MeV 

100 

hbsiS^Fs) 

762 ( 33772)7 10.706 3.555 

18.2 

1.88 X 10"^ 

10.853 

762 ( 23772)7 10.450 0.264 

2.0 

2.1 X 10“® 


762 ( 137 ) 2)7 10.165 0.053 

0.4 

4 X 10"3 


764(23 G4)7 10.770 

-2.781 

2.7 

2.8 X 10“® 


BB* 



33.2 MeV 

34.3 


B*B* 



58.2 MeV 

60.2 


BsBl 



3.32 MeV 

3.43 


BlBt 



1.93 MeV 

2.00 


Total 



96.7 MeV 

100 


TABLE XXII: Predicted partial widths and branching ra¬ 
tios for strong and electromagnetic transitions for the IG- 
wave states. For G-wave annihilation decays M designates 
^(»")( 0 )^ the fourth derivative of the radial wavefunction at 
the origin. See the caption to Table m for further details. 


Initial 

Final 

Mf 

M 

Width 

BR 

state 

state 

(GeV) 


(keV) 

(%) 

T5(1®G5) 

Xm(1®E4)7 

10.358 

3.057 23.1 

~ 100 

10.532 

T3(l®D3)7r7r 



2.2 X 10“^ 

9.5 X lO""^ 


Total 



23.1 

100 

T4(1®G4) 

Xi>4(l®F4)7 

10.358 

3.059 

1.4 

6.0 

10.531 

Xi>3(l®A3)7 

10.355 

3.032 

22.0 

94.0 


T3(l®D3)7r7r 



4.0 X 10“® 

1.7 X 10"^ 


T2(l®D2)7r7r 



1.8 X 10"^ 

7.7 X lO""^ 


Total 



23.4 

100 

T3(1®G3) 

Xi>4(l®F4)7 

10.358 

3.060 0.028 

0.12 

10.529 

Xi>3(l®T3)7 

10.355 

3.034 

1.8 

7.5 


Xf)2(l®F’2)7 

10.350 

3.005 

22.3 

92.4 


T3(l®D3)7r7r 



3.1 X 10"® 

1.3 X 10"® 


T2(l®D2)7r7r 



4.6 X 10“® 

1.9 X lO"'^ 


Ti(l®Di)7r7r 



1.7 X 10“"^ 

7.0 X lO"'^ 


Total 



24.1 

100 

77i,4(l^G4) 

hbsi^^ Fs)^ 

10.355 

3.034 

23.1 

~ 100 

10.530 

99 


1.005 

2.3 X 10~® 

1.0 X 10"^ 


?7i>2( 1^7)2 )7r7r 



2.2 X 10“"^ 

9.5 X lO"'^ 


Total 



23.1 

100 


best describes the properties of individual states but this 
underlines the importance of experimental input to test 
models and improve predictions. 

VII. SEARCH STRATEGIES 

An important motivation for this work is to suggest 
strategies to observe some of the missing bottomonium 
mesons. While there are similarities between searches at 
hadron colliders and e+e“ colliders there are important 
differences. As a consequence we will consider the two 
production channels separately. 


A. At the Large Hadron Collider 

1. Production 

An important ingredient needed in discussing searches 
for the missing bottomonium states at a hadron collider is 
an estimate of the production rate for the different states 
[SI ICTfSl] . The production cross sections for the T {nS) 
and XbJ states are in good agreement with predictions 
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TABLE XXIII: Predicted partial widths and branching ratios 
for electromagnetic transitions and OZI allowed strong decays 
for the 2G-wave states. Details of the OZI allowed decay 
amplitudes are described in the appendix. See the caption to 
Table HVl for further details. 


Initial 

Final M/ M 

Width 

BR 

state 

state (GeV) 

(keV) 

(%) 

T5(23G5) 

Xi,4(2®E4)7 10.622 3.598 20.6 

7.20 X 10“® 

10.772 

Xm(1®E4)7 10.358 0.186 

1.1 

3.8 X lO”'* 


BB 

25.9 MeV 

9.06 


BB* 

42.4 MeV 

14.8 


B*B* 

218 MeV 

76.2 


BsBs 

4.72 

1.65 X 10“® 


Total 

286 MeV 

100 

T4(23G4) 

Xi,4(2®E4)7 10.622 3.620 

1.3 

6.0 X lO""* 

10.770 

Xbsi'i^Fsh 10.619 3.570 19.7 

9.1 X 10“® 


BB* 

116 MeV 

53.7 


B*B* 

100. MeV 

46.3 


Total 

216 MeV 

100 

T3(23G3) 

Xi,4(2®E4)7 10.622 3.643 

0.025 

1.6 X 10"® 

10.769 

Xi>3(2®E3)7 10.619 3.593 

1.6 

1.0 X 10“® 


Xi>2(2®E2)7 10.615 3.538 

19.8 

1.28 X 10“^ 


BB 

10.3 MeV 

6.68 


BB* 

68.3 MeV 

44.3 


B*B* 

74.8 MeV 

48.5 


BsBs 

0.744 MeV 

0.482 


Total 

154.2 MeV 

100 

?7i,4(2^G4) 

hbai^^Fsh 10.619 3.573 20.7 

9.00 X 10“® 

10.770 

BB* 

108 MeV 

47.0 


B*B* 

122 MeV 

53.0 


Total 

230. MeV 

100 


of non-relativistic QCD (NRQCD) also referred to as the 
colour octet model. However we are interested in higher 
excitations with both higher principle quantum number 
and higher orbital angular momentum for which we are 
not aware of any existing calculations. To estimate pro¬ 
duction rates we use the NRQCD factorization approach 
to rescale measured event rates. In the NRQCD factor¬ 
ization approach the cross section goes like [5] 

a(if) = ^a„(A)(0^(A)) (II) 

n 

for quarkonium state H and where n denotes the colour, 
spin and angular momentum of the intermediate bb pair, 
(Jn{A) is the perturbative short distance (parton level) 
cross section and {O^ {A)) is the long distance matrix el¬ 
ement (LDME) which includes the colour octet QQ pair 
that evolves into quarkonium. We work with the assump¬ 
tion that the quarkonium state dependence resides pri¬ 
marily in the LDME which goes very roughly like (see 


Ref. [5T] 1 

(0(2^+iLy))cx(2J+I)^f® (12) 

where is the £ th derivative of the wavefunction at 

the origin and M is the mass of the state being produced. 
There are numerical factors, the operator coefficients of 
order I that have only been computed for the S- and P- 
wave states m- This gives, for example, an additional 
factor of 3 in the numerator for P-wave states in Eq. [T^ 
We note that at LO, NRQCD predictions are not in good 
agreement with experiment but at NLO the agreement 
is much better [5]. Some of the additional factors that 
contribute to the uncertainty in our crude estimates are 
not calculating the relative contributions of colour singlet 
and colour octet contributions, the neglect of higher order 
QCD corrections, the sensitivity of event rates to the pt 
cuts used in the analysis, and ignoring the dependence of 
detector efficiencies on photon energies. 

We will base our estimates on LHCb expectations but 
expect similar estimates for the collider experiments AT¬ 
LAS and CMS based on the measured event rates for 
bottomonium production by LHCb [31] , ATLAS |3S] and 
CMS |33|. However, there are differences between these 
experiments as LHCb covers the low px region while AT¬ 
LAS and CMS extend to higher px so that the production 
rates are not expected to be identical [321, only that the 
general trends are expected to be similar. 

To estimate production rates we start with the pro¬ 
duction rates measured by LHCb for LHC Run I and 
rescale them using Eq. Further, it is expected that 
the cross sections will more than double going from 8 TeV 
to 14 TeV [94] and the total integrated luminosity is ex¬ 
pected to be an order of magnitude larger for Run H com¬ 
pared to Run 1. LHCb observed ^ 1.07 x 10® T(IS')’s in 
the pAiJi~ final state for the combined 7 TeV and 8 TeV 
runs [34]. Taking into account the T(IS') —>■ p~ BR, 
roughly 4.3 x 10^ T(IS')’s were produced. Multiplying 
this value by the expected doubling of the production 
cross section going to the current 13 TeV centre-of-mass 
energy and the factor of 10 in integrated luminosity leads 
to 8.6 X 10® T(I5')’s as our starting point. We rescale this 
using Eq. [^ and use our estimates for the branching ra¬ 
tios for decay chains to estimate event rates. Our results 
can easily be rescaled to correct for the actual integrated 
luminosity. 

We find that this crude approach agrees with the LHCb 
event rates for nS and nP [34j within roughly a fac¬ 
tor of 2 , in some cases too small and in other cases too 
large. Considering the crudeness of these estimates and 
the many factors listed above which were not included 
we consider this to be acceptable agreement. We want 
to emphasize before proceeding that we only expect our 
estimates to be reliable as order of magnitude estimates 
but this should be sufficient to identify the most promis¬ 
ing channels to pursue. 

In the following subsections we generally focus on bb 
states below BB threshold, as the BR’s for decay chains 












20 


TABLE XXIV: Summary of lowest order expressions and first order QCD corrections with Os computed at the mass scale of 
the decaying state (see Section IV for references). 


Process 
^So gg 

^Sq 77 
->■ ggg 
'Si ^ 751? 

^Si ->■ 777 

->• e+e- 
'-P 2 -S' gg 


Rate 




~^Q 
40(7r^-9)a5 


|vl/(0)| = 


32(71^ — 9 )eQCiCi^ 


9m- 


Q 


16(7r^-9)e^a 


Zmj. 


|d/( 0 )| = 
-|d/( 0 )p 


Q 
2 2 


'P 2 —>■ 77 
^ gg + g 
'-Po -S' gg 




^^lR'„p(0)l^ln(mQ(R}) 


^IKp( 0 )| = 


'fo -t 77 

'P ’1 ^ ggg 

'P>3 -S' ggg 
'P>2 -t 555 
'P>1 -t 555 
e+e- 

'P >2 -S' 55 

'P 4 -S' 55 
'p’s -t 55 
'P ’2 -S' 55 

'<^4 —S' 55 


^^IKp(o)|' 

i|Kp(0)pln(mQ(R)) 


97rm® 

lOck^ 

97rm® 

760a 


|R"o(0)pln(4mQ{R)) 

|R::i,(0)pln(4mQ{R)) 


^|i?:(i,(0)|^ln(4mQ{R)) 




mb 

2al 


27m^ 


27m» 

919^1 

135m^ 


3-7rm^ 


'nD( 0 )|" 

Kd(o)|' 

|i?"p( 0 )|' 

Kp( 0 )|' 

Kp(0)|' 

IKg(o)I' 


Correction 

(1 + (for bl) 

( 1 -^) 

(1 - (for bb) 

(1 - P^) (for 6 &) 

(■\ _ 12 . 60:3 \ 


(1-^) (X<>2(1P)) 

(1 + ^) (Xs>2(2P)) 
(1 + (Xs>2(3P)) 

(1 + ^) (Xs>2(4P)) 
( 1 -^) 


(1 + ^) (Xw(lP)) 
(1 + i^^) (Xs>o(2P)) 
(1 + (Xm(3P)) 

(1 + i^) (Xs>o(4P)) 

(1 + P^) 


( 1 -^) 


originating from states above BB threshold will generally 
result in too few events to be observable. Likewise, the 
production cross sections for high L states are suppressed 
by large powers of masses in the denominator. 

We also focus on decay chains involving radiative tran¬ 
sitions although hadronic transitions with charged pi- 
ons often have higher detection efficiencies. However, 
hadronic transitions are not nearly as well understood as 
radiative transitions so we were only able to even at¬ 
tempt to estimate a limited number of cases that, as 
discussed, we could relate to measured transitions. In 
addition, in many cases of interest, hadronic transitions 


are expected to be small. Nevertheless, as demonstrated 
by the study of the T{l^Dj) by the BaBar collaboration 
[H] and the hb{lP) and hb{2P) by the Belle collaboration 
m, hadronic transitions offer another means to hnd and 
study bottomonium states. We include a few examples 
in the tables that follow but they are by no means an ex¬ 
haustive compilation and we encourage experimentalists 
to not neglect this decay mode. 
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2. The 3S and higher excited S-wave states 

We start with the S'-wave states. Our interest is that 
they will be produced in large quantities and their decay 
chains include states we are interested in such as excited 
P- and H-waves and will therefore add to the statistics 
for those states. We therefore focus on decay chains that 
include these states. 

The 3S decay chains of interest are given in Ta¬ 
ble IXXVIl The estimates for the number of events ex¬ 
pected for LHCb are included but we also include a col¬ 
umn with estimates for e"'"e“ collisions expected by Belle 
II which we discuss in Sec. Ivira What is relevant is 
that numerous \D states will be produced in this man¬ 
ner and when added to those produced directly and via 
P-wave initial states will give rise to significant statistics 
in 27 -I- final states. In addition, it might be pos¬ 
sible to observe the r]{,{3S) in a final state from 

the %(3S') —>■ T(I5') MI transition. 

The 4^ S'! state can decay to 3P states which can sub- 
seq uently d ecay to 2D or ID states, as shown in Ta¬ 
ble XXVII The is above the BB threshold so has 
a much larger total width than the lower mass 5'-waves 
leading to a much smaller BR for radiative transitions. 
Decay chains to P-waves involve too many transitions 
making them difficult to reconstruct so we do not include 
them in our tables. We also include decay chains which 
might be of interest for e+e“ studies but would result 
in insufficient statistics to be relevant to hadron collider 
studies. We do not include the 4^S'o state as the decay 
chains have too small combined BR to be observed. 

For the the BR’s to the 4^Pj states are 0(10“^) 
and the BR’s of the 4^Pj are ^ 10“"^ so this product 
BR is quite small. When we include BR’s to interesting 
states such as the 3^Dj states the product BR’s are likely 
to be far too small to be observable. The BR’s to the 
3^Pj states are comparable to the 45 —)■ 3P transitions, 
0(10“^) so it might be possible to see 3P states start¬ 
ing from the It is not likely that the 2D and ID 

states can be observed in decay chains originating from 
the 5^S'!. We arrive at similar conclusions for the 6^Si 
and conclude that the only possible states that might be 
observed are the 3^Pj states. 


3. The 2P states 

The 2P states are of course well known. We include 
them as they can decay to the ID states and hence con¬ 
tribute to the ID event rates which was the discovery 
channel for the I^Zl 2 state [SS] (see also Ref. [M]). We 
only include event chains relevant to these final states 
which we list in Table |XXVIII[ In an attempt to re¬ 
duce the theoretical uncertainties in the 2P ID 
BR’s, rather than using the predictions for total widths 
and BR’s in Table IX we estimate the Xf,j(2P) total 
widths using the PDG values |31] for the BR’s for the 
2P 2S + ^ and 2P —>■ IP -|- 7 transitions together 


with our predictions for these partial widths as was 
described in Section VB. The values for r( 2 ^P 2 ) and 
r(2^Pi) were given there. Similarly we find r(2^Po) = 
247 ± 93 keV. Combining these total widths with the 
par tial widths for 2P —>■ ID transitions given in Ta- 
r( 23 p 2 ^ I^D 3 . 2 ,i) = 1.5, 0.3 and 0.03 keV 


IX 


ble 

respectively, we obtain the corresponding BR’s of 1.2%, 
0.2% and 0.02%. Likewise r{2^Pi l^D 2 ,i) = 1.2 and 
0.5 keV give the corresponding BR’s 1.9% and 0.80%, 
and r(23Po ^ l^Di) = I.O keV has BR = 0.4%. There 
is uncertainty in these estimates as can be seen by com¬ 
paring the PDG values [32] we used in our estimates to 
recent BaBar measurements nnni and by comparing our 
predictions for the partial widths to those of Kwong and 
Rosner [3S]. Nevertheless these estimates are sufficient 
for estimating the 2P —>■ ID + 7 BR’s and the result¬ 
ing event rates for the purposes of identifying promising 
channels. The important point is that in all cases, signifi¬ 
cant numbers of bottomonium D-waves will be produced 
from P-wave production and decay which will improve 
the statistics from those that are produced directly. 


4-. The 3P states 

The observation of the Xb{3P) by the ATLAS col¬ 
laboration [T] through its radiative transitions to T(IP) 
and T(2P) with T(IP, 2S) -T was the first new 

particle discovered at the LHC. This decay chain rep¬ 
resents a clean experimental signature with the two fi¬ 
nal state muons a clean signal to trigger on. The 
Xb{3P) was confirmed by the DO collaboration |I()Ij 
and by the LHCb collaboration |33l [34]. Further, 
LHCb identified the state as the Xbi{3P) with mass 
I05I5.7^3;g(stat)t2'i(syst) MeV/c^ [33] . Using the ap¬ 
proach outlined above we calculate for Run I ^ 243 
events for the decay chain x(3P) —>■ T(I5') —>■ 

~ 371 for the T(2S') decay chain, and ~ 1030 for the 
T (35') decay chain compared to the observed numbers of 
events of 329±59, I2I±3I and I82±23 respectively [34] . 
The agreement is reasonable for the T(I5) chain but be¬ 
comes decreasingly so for the T(25) and T(35) chains. 
We assume this is due to decreasing photon detection 
efficiency as the photon energy goes down. Overall the 
agreement is not unreasonable considering the approxi¬ 
mations used to obtain these values and our neglect of 
detector efficiencies and helps the reader judge the gen¬ 
eral reliability of our predictions for event rates. 


In Table XXIX we summarize the event rates expected 


for the most promising decay chains. The Xboi^P) state is 
much broader and decays predominantly to light hadrons 
via gluon intermediate states with a small BR to 7 T(I 5 ) 
(^ 10 “"^) so it would be quite challenging to observe al¬ 
though it might be possible to observe via hadronic tran¬ 
sitions to the I^P 2 state which subsequently undergoes a 
radiative transition to the I^5i state. 

For both the Xb2{3P) and Xbi{3P) the decay chains 
to T(n5) —>■ where n = I, 2, 3 give rise to the 
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largest event rates and are the simplest to reconstruct so 
it is not surprising that these were the discovery chan¬ 
nels. Run II should provide sufficient statistics to sep¬ 
arately fit the Xb 2 {‘iP) and Xbii^P) using these decay 
chains. Other decay chains are potentially interesting as 
they involve the undiscovered 2 ^D 3 , 2 ^D 2 , 2^Di, l^Ds, 
1^F4 and bb states. However they generally have 

multiple photons in the final state making it difficult to 
reconstruct the initial particle. In addition, some of the 
photons in the decay chains are relatively low energy so 
could have low detection efficiencies. The one exception 
is 3^Pi 2^Dix —>■ IJF ■ This final state is relatively 
clean but has a low combined BR. As we noted, we would 
not be surprised if our production rates are off by an or¬ 
der of magnitude so we do not rule out the possibility 
that the 2^Di could be observed in this process. 

We do not include decay chains for the /it,(3P) state, 
with two exceptions, as we believe it would be very dif¬ 
ficult to reconstruct in a hadronic environment for the 
following reasons. It decays predominantly to hadronic 
final states 70% of the time) which would be diffi¬ 
cult to identify in hadron collisions. The next largest 
BR’s are to the r]b{3S), ribi2S) and r]b{lS). All of these 
decay almost 100% of the time to hadronic final states. 
The r]b{3S) and r]b{2S) have small BR’s to hb{nP) and 
very small BR’s to T(15'). While there are several hun¬ 
dred hbS produced they would be difficult to see and 
the combined BR to T(15') —>■ iJp fJT is too small to pro¬ 
duce sufficient numbers to observe. We include the decay 
chain hb{3P) r? 6 (lS ')7 —>■ gg as it might be possible to 
reconstruct the hb{3P) using the one final state photon 
and possibly simple r]b hadronic final states and likewise 
the decay chain hb{3P) —>■ hb{lP) -I- tttt —>• r]b{lS)x —>■ gg 
for the same reasons. 


5. The 4P and 5P states 

The 4P and 5P states are above BB threshold and 
have total widths of the order of tens of MeV so that BR’s 
for radiative transitions are relatively small, 0(10“"^). 
However they undergo radiative transitions to T(nS') 
states which decay to jjr offering a clean final state to 
study. Our estimates for the expected number of events 
from these decay chains are given in Table [XXX| While 
the number of expected events is small it is possible that 
our estimates are off by an order of magnitude. Also, our 
estimates are based on using LHCb event numbers and 
ATLAS and CMS with their different capabilities might 
be able to observe more events. Thus, for completeness 
we include estimates for both the AP and 5P states. The 
most important message is that the higher energy and lu¬ 
minosity of LHC Run B could potentially observe some 
of the 4P and 5P states which would be an important 
test of models and Lattice QCD results. 


6. The ID states 


The production rate for the P-waves is significantly 
lower than that of P-waves. This is a consequence of 
the LDMEs going like |p(')( 0 )| VM2'+2 so the produc¬ 
tion cross section is significantly suppressed by the mass 
factor in the denominator. The mass factor is simply a 
consequence of the dimensionality of the Zth derivative of 
the wavefunction and will increasingly suppress the cross 
sections going to higher I multiplets. 

Nevertheless, it is expected that the ID states can be 
produced in sufficient quantity to be observed. There are 
three sources of the D states; direct production and from 
decay chains originating with the 3^Si and 2^Pj states. 
The decay chains and estimated number of events from 
direct production are given in Table |XXXI[ We estimate 
that direct production of the P-waves will yield ~ 100 
T 3 ( 1 P), ~ 150 T 2 ( 1 P) and ~ 50 Ti(lP) events. This 
is roughly comparable to the number of 3P events ob¬ 
served by LHCb. There are two differences. The 3P 
events were comprised of one photon and a g~ pair 
while the IP events are generally comprised of two pho¬ 
tons and a g'^ g~ pair making them more difficult to 
reconstruct. On the other hand, we can estimate the 
photon energies fairly accurately because the T 2 ( 1 P) 
mass has been measured. An exception is the Ti(lP) 
which can decay directly to a g'^g~ final state. Un¬ 
fortunately the BR appears to be too small to produce 
a sufficient number of events to find this state in this 
channel. On the other hand, as we have pointed out a 
number of times, our estimates can easily be off by an 
order of magnitude. In addition to direct production the 
IP states will also be produced via transitions originat¬ 
ing with 3^Si and 2^Pj. In fact, cascades originating 
from the 3^Si contribute the largest number of events to 
the 2xg'^g~ signal with roughly another 20 % originat¬ 
ing from 2^Pj production. Specifically we expect ^ 2700 
events in l^Pa — 7 l^P 2 — 77 l^ 5 'i —>■ 77 /^+^“, ^ 6500 
in 1 ^P 2 —>■ 7 l^Pi —?► 77 l^S'i —>■ ~ 1200 

in l^Pi —>• 7 l^Pi —)■ 77 l^S'i —>■ jxg'^g~. Other decay 
chains will contribute to l^Pj production but they will 
have different photon energies so we only give estimates 
for the decay chains with the largest statistics. 

We have also included the decay chains l^Pj —>■ 
tt+tt” —>• g'^g~ and 1 ^P 2 —>■ l^S'o + 7r''‘7r“ —>• gg despite 
yielding few events as they offer a signal complementary 
to those involving photons. 

Finally, we mention that the rib 2 will decay predomi¬ 
nantly via rib 2 {^D) jhb 77776 ( 15 ') with the ijbi^S) 
decaying to hadrons. We do not see how the 77 ^ 2 (IP) 
can be reconstructed given the hadronic final state but 
perhaps experimentalists will come up with a clever ap¬ 
proach that we have not considered. 
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7. The 2D and higher D-wave states 


For the 2D states in Table XXXII we focus on the 
simplest decay chains with 2 photons and a 
the final state which are most likely to be reconstructed. 
To estimate the number of expected events, we include 
contributions to the 77 / 1 +^“ final state from direct pro¬ 
duction and from decay chains originating from 3^Pj pro¬ 
duction. There are a number of possible decay chains but 
the ones with the largest expected statistics are ~ 125 
events for 2 ^ 1)3 —)• ^2^P2 —t 772 ^S'i —)■ 77 ^+/^“, ~ 252 
events for 2 ^D 2 —t 72 ^Pi —>■ ^^2^Si -T 77 /^+^“ and 
~ 37 events for 2^I?i — 'y2^Pi —>• 772 ^ 6 '! — 77 /^+^“. 
In all cases more events are expected from 3^Pj produc¬ 
tion than from direct 2D production. 

For the 3D and AD states we expect that each member 
of both multiplets will have of the order of 10 "^ produced. 
However the predicted widths are 0(100 MeV) so the 
BR’s for radiative transitions will be small. Thus, we 
only expect that they can be observed in BB, BB* or 
B*B* final states if they can be reconstructed with high 
enough efficiencies and separated from backgrounds. 


8. The nF states 

The production rate decreases quite dramatically for 
states with larger L as a consequence of our estimates 
which use Eq. where the NR approximation of the 
cross section goes like the /th derivative of the wavefunc- 
tion at the origin with the corresponding mass in the 
denominator needed for dimensional reasons. We only 
expect that ~ 200 for each of the IF states will be pro¬ 
duced. Once the BR’s for the decay chains are included 
we expect that only 1 or less events will result. Con¬ 
sidering experimental challenges in making these mea¬ 
surements we do not expect that the IF states will be 
observed from direct production. Nevertheless, we in¬ 
clude the dominant decay chains in Table |XXXIII| for 
completeness. 

We do expect a small number of IF states; ^ 28 l^F 4 ’s, 
~ 44 I^Fa’s and ^ 2 1 ^F 2 , to be produced via radiative 
transitions originating with 3^Pj and 2^Dj states. Since 
there are 3 7 ’s in the final state it is unlikely that the IF 
states will be observed in hadron production. 

We expect the 2F and 3F multiplets to be even more 
challenging to observe using radiative transition decay 
chains primarily because they are above BB threshold 
and are therefore broader, ranging from 2.8 MeV for 
the 2 ^F 4 state to 88.6 MeV for the 2 ^F 2 state. These 
states are very close to BB and BB* threshold and there¬ 
fore are very sensitive to available phase space. If our 
mass predictions are too high it is possible that the to¬ 
tal widths could be significantly smaller leading to sig¬ 
nificantly larger BR’s to the decay chains we have been 
focusing on. Nevertheless, given our expectations for the 
IF states we do not feel it is likely that they would be 
discovered using radiative transitions. If the B and B* 


mesons can be observed with high efficiencies the excited 
F-wave states might be observed in BB and BB* final 
states. 

The 3F multiplets are sufficiently above BB and BB* 
threshold that they are much broader. The only possi¬ 
bility that they might be observed would be in BB, BB* 
or B*B* final states. 


9. The nG states 

Even more so than the F-waves, the G-wave produc¬ 
tion cross sections are highly suppressed so we expect 
only 0{1) meson will be produced for each of the states 
in the G-wave multiplets in Run II which is far too small 
a number to have any hope of being seen. 


B. At e''"e Colliders 


The question we wish to address in this section is 
whether previously unobserved states can be observed in 
collisions and outline how to do so. We will focus 
on 1 states, the n^Si and n^Di states, since only 1 
states can be produced directly in e“''e“ collisions. To 
estimate the number of events requires a detailed Monte 
Carlo study that includes beam spread and initial state 
radiation which is beyond the scope of this work. Instead 
we will use cross sections derived from Belle and BaBar 
events along with integrated luminosities suggested for 
Belle II. 

We start with the T(35') where we use the cross sec¬ 
tion of 4 nb based on Belle and BaBar measurements 
|I 02 j . An integrated luminosity of 250 fb“^ would yield 
10® T(3S')’s, about a factor of 7 1/2 times the combined 
Belle-BaBar dataset. We give estimates of number of 
events in Table IXXVIl based on this but the numbers can 
be easily rescaled for other integrated luminosities. We 
expect that G( 10 ®) l^Fj’s will be produced via radiative 
transitions for each J value, which should be sufficient to 
identify each of the l^Dj states and measure and com¬ 
pare their masses to theoretical predictions. 

The expectations are to accumulate several ab“^ of 
integrated luminosity at the T(4S'). The cross section 
of ^ 1 nb would yield several billion T(4S')’s depend¬ 
ing on the eventual integrated luminosity. We estimate 
the number of events given in Table |XXVII assuming 
that 10^° T(45')’s will be produced based on 10 ab“^ 
of integrated Inminosity but these numbers can easily 
be rescaled assuming different values of integrated lumi¬ 
nosities. We expect sufficient events in the chains that 
proceed via the 3 ®F 2 and 3®Fi so that these states should 
be observed in radiative decays of the T(4S'). We don’t 
expect that the 3®Fo will be observed in this manner. 
Another interesting possibility for studying the 3P states 
via radiative transitions from the T (dF) utilizes hadronic 
transitions from the T(35') and T{2S) to T(2S') or T(IS') 
in the decay chain. This is experimentally very clean and 
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would yield some tens of events for 3 ^P 2 and 3^Pi inter¬ 
mediate states but only 0(1) event for the 3^Po- This 
might be sufficient to resolve some of the J states. It 
would have been interesting to be able to observe the 2D 
states in radiative transitions originating with the T (4S') 
but this does not appear to be likely. Of all the 2^Dj 
states the decay chain proceeding via the 2^Di state will 
have the largest statistics although not sufficient to be 
observed considering the 4 photons in the decay chain. 

It is unlikely that sitting on the T(5S') or T( 6 S') would 
produce enough bottomonium states via radiative transi¬ 
tions to be seen. This stems from the T(5S') and T( 6 S') 
having large total widths which leads to small BR’s to 
e“'"e“ and subsequently smaller production cross sections 
than the lower mass ^Si states. This also results in small 
BR’s for photon decays. For the T{5S) the Belle collab¬ 
oration measures its e“'"e“ production cross section to be 
0.3 nb |102j . The radiative transition with the largest 
BR is to 4 ^P 2 with BR = 1.6 x 10~^ y ielding ~ 5 x 10^ 
4 ^P 2 ’s per 100 fb“^. From Table XXXj the combined BR 
for representative 4 ^P 2 decay chains are 0 ( 10 “®) so it is 
unlikely that bb states could be seen in decays originating 
from the T(5iS') in e+e” collisions. We repeat the exer¬ 
cise for the T( 6 S'). We assume an e+e“ production cross 
section similar to that of the T{5S) although in fact it 
should be smaller due to its smaller BR to e+e“. For the 
T( 6 S') the largest BR is expected to be to the 3®Pj states, 
0(10“®). However, in this case the combined BR’s for 
some interesting 3P decay chains are 0(10“"*). Putting 
this together we expect ^ 0.5 event/100 fb“^ via the 
3 ®P 2 intermediate state. While still below an event rate 
needed to study these states it does offer some hope given 
the uncertainties in our assumption for the production 
cross section and with higher statistics. 

Observation of the n^Di states in e'''e“ collisions is 
interesting, even more so if they can be produced in suf¬ 
ficient numbers to see previously unobserved excited bot¬ 
tomonium states in their decays. Unfortunately the pro¬ 
duction cross section is proportional to the BR to e“'"e“ 
which is roughly three orders of magnitude lower than 
for the S'-wave states. The small number of signal events 
will also make it challenging to see the n^Di states over 
backgrounds. Given these caveats we make a rough es¬ 
timate of the number of rt’Di produced by multiplying 
the ratio of nD/2S BR’s to e+e“ times the 2^51 cross 
section to obtain the n^Di production cross section. For 
the l^Di this gives ~ 13 pb. For 100 fb“^ of integrated 
luminosity we expect ~ 1.3 x 10® l^Pi’s to be produced. 
The expected number of events for several l^Di decay 
chains is given in Table pQCXI| This will yield 50 events in 
l^Di —>■ but many more in radiative decay chains 

via intermediate IP states. 


In the same way as we estimated CT(e“''e“ —?► l®Pi) we 
obtain a{e~^e~ 2®Pi) ^ 18 pb. For 100 fb“^ of inte¬ 

grated luminosity we expect ^ 2 x 10 ® 2 ®Pi’s produced. 
We use this number to estimate the number of events in 
the different decay chains given in Table |XXXIV[ With 
these assumptions it should be possible to observe the 


2®Pi. While the mode may be the cleanest, other 

modes proceeding via radiative transition decay chains 
offer higher statistics, in particular decay chains proceed¬ 
ing by IP and 2P intermediate states. More interesting 
is the possibility of observing the 1 ®P 2 for the hrst time. 
The decay chain with the highest statistics proceeds via 
1 ®P 2 —>■ I®Pi —)■ l®Pi which would also yield informa¬ 
tion on the l®Pi state. One can easily scale our projected 
events up or down to reflect actual cross sections and in¬ 
tegrated luminosities. 

The 3D and AD are above BB threshold so have total 
widths 0(100 MeV) resulting in BR ^ 10“® and a cross 
section ~ 8 fb. For reasonable integrated luminosities 
and u sin g the B R’s for the 3®P>i and 4®P>i given in Tables 
XVI and XVIII we do not expect that the 3®Pi and 4®Pi 
states can be observed in e+e“ collisions, as suggested by 
the results given in Tables XXXV and XXXVI| 


VIII. SUMMARY 

In this paper we calculated the properties of bottomo¬ 
nium mesons including masses, radiative transitions, an¬ 
nihilation decays, hadronic transitions and strong OZI 
allowed decays for states above threshold. These re¬ 
sults were included in extensive tables with estimated 
BR’s to different Hnal states. While we are interested in 
how these predictions fare against experimental measure¬ 
ments as a test of our understanding in the context of the 
constituent quark model, the main objective of this work 
is to make predictions that can assist experimentalists in 
finding missing bottomonium states and measuring their 
properties. 

We estimated the number of events expected in Run 
II of the LHC in the context of the LHCb experiment 
but they should also be relevant to the ATLAS and 
CMS experiments. We expect that significant numbers 
of Xb2{3P) and Xbi(3P) will be produced and decay via 
radiative transitions to the T(IS') and T(2S') which will 
subsequently decay to ■ Likewise we expect a large 
number of X62(4P) and Xbi(4P)’s to decay to T(3S'). 
Thus, a promising search strategy for excited P-wave 
mesons is to reconstruct T’s in and look at the 

invariant mass distributions of the T’s with one 7 . We 
also expect that r]i,{3S) can be seen in this final state. 

Turning to the P-waves, the l^D^ and 1 ®P 2 will un¬ 
dergo radiative transitions to the 1 ®P 2 and l®Pi respec¬ 
tively so they might be seen in hnal states with 77 /i+/i“. 
Similarly, the 2D states will decay to 2 ®P 2 _i which decay 
to 2^Si- Thus, it might be possible to see most of the ID 
and 2D spin triplet multiplets in the 77 /i’^/r“ hnal state. 
A challenge is that the photon energy for 2 ®P 3 —>■ 2 ®P 2 
is almost identical to that from the 2 ®Pi — 2 ®Pi tran¬ 
sition so that one would need to be careful in looking at 
the invariant mass distributions of the hnal state. 

For e“'"e“ collisions, large numbers of l®Pj states will 
be produced by sitting on the T {3S) so it might be pos¬ 
sible to resolve the three states and determine the split- 

















25 


tings between members of the multiplets. Sitting on the 
T(4S') will produce 3 ^P 2 ,i in radiative transitions. It 
should be possible to observe the l^Di and 2^Di by an 
energy scan at the appropriate energy. Sitting on the 
2^Di resonance it might be possible to observe the 1 ^F 2 
via radiative transitions from the 2^Di. 

The LHC experiments and Belle II hold the promise 
to increase our knowledge of bottomonium mesons. This 
improved knowledge will test the reliability of models of 
quarkonium physics. Lattice QCD is making ever more 
precise calculations of bottomonium mesons and it is im¬ 
portant that these calculations be held to account by 
experiment. We expect that the phenomenological pre¬ 
dictions presented in this paper will be a useful tool for 
experimentalists to do so. 
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Appendix A: The ®Po Model 

The quark pair creation model [niisiiTiiniiss] 
is used to calculate OZI allowed strong bottomonium 
decays. Given that details of the calculations such as 
phase conventions are important and not always clearly 
stated in the literature we summarize the details of the 
model to assist an interested reader in reproducing 
our results. 

In the quark pair creation model [niiiiiaiiTi 
I85j . a qq pair is created from the vacuum in a ^Pq 
state (the quantum numbers of the vacuum). The an¬ 
gular momentum and spin of the created qq pair are 


therefore Lp = Sp = 1, and Jp = 0, so that 
Mpp = —Msp = m. The transition operator for qq 
pair creation can be written as 

T = — 37 ^( 11 ; m — to|00)-\/9(^ J d^pq(fipqd^{pq + Pq) 

where b^q{pq) and dq{pq) are the creation operators for the 
quark and antiquark, respectively. The momenta of the 
created quark, pg, and the created antiquark, pq, are in¬ 
tegrated over all possible values, such that the delta func¬ 
tion ensures that their total momentum is zero in their 
centre-of-mass frame. The spin triplet state of the cre¬ 
ated qq pair is described by its spin wavefunction Xi-m 
and the momentum-space distribution of the created pair 
is described by the solid harmonic, written in terms of the 

spherical harmonic as (^) = ( 6 *fc,'/'fc)- The 

SU{3) flavour singlet wavefunction of the created pair is 
(j)o = -^{uu+dd+ss) and its colour singlet wavefunction 
is loq. The overall factor of 3 in Eq. |A1| will cancel out 
when evaluating the colour overlap. The factor of •\/ 967 r 
arises from the normalization and field theory conven¬ 
tions of Refs. [II1IH3IH7]. The amplitude for quark-pair 
creation from the vacuum can therefore be described by a 
single free parameter, 7 . Any differences in the constant 
factors that appear in Eq. |Al| simply result in a rescaling 
of the value of 7 . For example, the value of 7 used in 
Refs, [m [89] is larger than ours by a factor of V 967 r due 
to the absence of this factor in their T operator. 

The S-matrix for the meson strong decay A —>■ BC is 
defined as 

S = I- 2-Ki5{Ef - E,)T (A2) 

so that 

{BC\T\A)= S^pA - PB - Pc)M^-'^’^-’b^-’c (A3) 

where, using the normalization from Refs. [nilHSlIHT], 
the helicity amplitude is given by 


= 7 ^ {LaMl^SaMs^\JaMjP{LbMlpSbMsp\JbMjp){LcMl^ScMsc\JcAIji^) ( A 4 ) 

(1ml - rn\00){xs^pMsgXTcMsc IXs^Ms^Xi-m) 

('/'bVc <(>o'‘)-f(4wi,m2,m3) -b (-1)^+'®^+'^^+'®'^((()f (()^'‘|(/)^2(/)^'^)/(-P,m2,mi,m3) 


where the sum is over Ms^, Mp^, Ms^, 

Msc , and m. The two terms in the last factor correspond 
to the two possible diagrams. In the first diagram, quark 


1 from meson A ends up in meson B and antiquark 2 from 
meson A ends up in meson C. In the second diagram, 
quark 1 from meson A ends up in meson C and antiquark 
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2 from meson A ends up in meson B. Indices 3 and 4 The momentum space integral for the first diagram is 
refer to the created quark and antiquark, respectively. given by 


/(P,mi,TO 2 ,m 3 ) = J d^pyim{p)tpUA,LA,ML^ {p + p) '4 ’ub,Lb,Mlj, (p + 


I* 

\^+ 

m3 

P] 


riB 

mi + m3 

) 


;* 

(p + 

m3 


(A5) 

nc i-L'C 

m2 - 1 - m3 


where mi, m 2 and m 3 = 1114 are the constituent quark 
masses and we have taken P = Pb = —Pc in the centre- 
of-mass frame oi A. To evaluate the spatial integral 
for the second diagram, we simply interchange B -tCf C^ 
which amounts to m aking the replacement mi -f-)- m 2 and 
P —>■ — P in Eq. A5 leading to the second term of Eq. 


A4 


The techniques found in Appendix A of [lOSj were use¬ 
ful in simplifying and evaluating the spatial integrals. Eor 
the meson space wavefunctions, we use the momentum- 
space simple harmonic oscillator (SHO) wavefunctions, 
given by 


V'nLMi {P) = PnL {p)YlMb [dp, (j^p) 
where the radial wavefunctions are given by 


(A 6 ) 


R 


SHO 


(P) = 


{-in-tY 


2 n 


r(7z L 2 ) 




(A7) 


The colour matrix element 




(AS) 


does not explicitly appear in Eq. Ad since it cancels the 


overall factor of 3 from Eq. A1 


The flavour matrix element can be easily found by writ¬ 
ing the flavour wavefunctions of mesons A, B and C and 
that of the created quark pair as 5 x 5 matrices with rows 
indicating the quark flavour (u, d, s, c, b) and columns in¬ 
dicating the anti-quark flavour (u, d, s, c, b). For example, 
the flavour wavefunction for the created qq pair is 


/ 1 0 0 0 0 \ 

^ 0 10 0 0 
dd + ss) = ^ 0 0 10 0 

0 0 0 0 0 
Voodoo/ 

The flavour overlaps for each of the two terms in Eq. 
are therefore given by 




(A9) 


Z/+ — 

and Ln ^ (p^//3^) is an associated Laguerre polynomial. 
We use the SHO wavefunctions such that a meson with 
quantum numbers in spectroscopic notation 

uses ' 4 ’Y.YYlMb momentum-space wavefunction. 

The values we use for the effective harmonic oscillator 
parameter, /3, are listed in Tables m 


= Pr[(l)l(j}B4>0 (l^c] (AlO) 

= T'r[(j)\(l)c4'o (I^b] (All) 

The spin matrix elements for the first and second di¬ 
agrams are written in terms of the Wigner 9j symbols 
PM 1 as 


(XSaMsaXScMg^ IXSAMs^Xl-m) — ( —l)^'*"'^'^ y/3{2SA + 1)(2S'b + l)(2S'c -|- 1) 


5 5 


{SbMsbScMsc\SMs){SaMsa^ - m\SMs) < 

1 

2 

\ 1 

\ (A12) 

S,Ms I 

Sb 

Sc S 

J 

/•v32 34 \ _ / i U-I-Sa-I-Sb+Sc 32 i 12 34 \ 

[XSBMsgXScMs^lXSAMs^Xl-ml “ I.--*-! USsMs^ XScMg^, IXs^Mg^ Xl-ml 



(A13) 


where the spin matrix element for the second diagram 
was obtained using an alternative definition for the 9j 
symbols that couple the quarks differently |104) . This 


expression, given in Eq. A13 was used to simplify Eq.|A4| 


Using the Jacob-Wick formula |105l 1106 
amplitudes ^ given by Eq. 


the helicity 


A4 


are con- 
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verted to partial wave amplitudes via 


M^^{P) = 


^Att{2L + 1 ) 
2Ja + 1 


E 

^ Q 


{LOSMj ^I JaMj^){JbMj, JcMj^\SMj^)M^^^^(Pz) 


Mja —^Jb '^^Jc 

(A14) 


where S = Jb + Jc and Ja = L + S such that 

\Jb — Jc\ < S < Jb + Jc (A15) 

\Ja-S\<L<Ja + S (A16) 

and the outgoing momentum of meson B, P = Pz, is 
chosen to lie along the z-axis in the centre-of-mass frame 
of meson A so that the helicities and angular momentum 
projections are related by Mj^ = A^, Mj^ = —Xq and 
Mj^ = Mj^ +Mj^ = Xb — Xc- This on-shell momentum 
is conveniently written in terms of the masses of mesons 
A, B and C as 


where Eb = a/M^ -|- P'^, Eq = \/Mq + and 5 is a 
symmetry factor given by 


5 = 


1 

1 + Sbc 


A if i? and C are identical 
1 otherwise 


(A19) 


Finally, the strong decay width for a given decay mode 
of meson A is just the sum of its partial widths: 


- {Mb + Mc)2] [Ml - (Mb - M^] 

2Ma 

(A17) 

Using relativistic phase space, as described in Ref. [HI 
I85| , the partial width for a given partial wave amplitude 
is given by 

= 2TrPS ^^^^ (A18) 

Ma 


r = ^r^^. (A20) 

L,S 


The calculations of the strong decay widths, as outlined 
in this section, were performed using the Mathematica 
software package, version 7.0 m- 
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TABLE XXV: Expressions for hadronic transitions in terms 
of the reduced amplitudes Ak{i' ,£). Note that reduced ampli¬ 
tudes are dependent on the initial and final states. Because 
the Ai contributions are suppressed we follow the usual prac¬ 
tice and will take Ai{i',i) — 0 although we include them in 
the table. The details on how we obtain numerical estimates 
for amplitudes are described in the text. 


Process 



Expression 



Si -I- TTTT 


Ao( 0 , 0 ) 


'S'o 

So -b TTTr 


Ao( 0 , 0 ) 


3 p 2 ^3 

P2 -b TTTT 

1 Ao(l, 1 ) -b 

2 Ai( 1 , 1 ) -b 

^^2(1,1) 

3 p 2 ^3 

Pi -b TTTT 

3 ^Ai( 1 , 1 ) + |,A 2 ( 1 , 1 ) 

3 p 2 _^3 

Po -b TTTT 


11^2(1,1) 



P2 -b TTTT 

36 ^ 1 ' 

( 1 , 1 ) + |A 2 ( 1 , 1 ) 


Pi -b TTTT 

|Ao( 1 , 1 ) -b 

1 ) -b 

1 ^ 412 ( 1 , 1 ) 

3 Pi 

Po -b TTTT 




"Po 

P2 -b TTTT 


1^2(1, 1 ) 


"Po 

Pi -b TTTT 




"Po 

Po -b TTTT 




IPl 

Pi -b TTTT 

|Ao(l, 1 ) 4 

- |Ai( 1 , 1 ) -b 

5412(1,1) 


P3 -b TTTT 

Mo( 2 , 2 ) + 

r5M2,2) + 

^^2(2,2) 


P2 + TTTr 

^Ai(2,2) +^^2(2,2) 


Pi -b TTTT 


145^2(2,2) 



P3 -b TTTT 

225^1' 

(2,2) +^^2(2,2) 
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i^i(2,2) + 
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Pi -b TTTT 

100^1 (2’ 2) 100^2 
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P3 -b TTTT 
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P2 -b TTTT 
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|bl 2 ( 3 ,l) 
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TABLE XXVI: The 3S Decay chains, branching ratios and event estimates for LHCb Run II and Belle II. The number of events 
in the pp column are based on producing 3.1 x 10® T(3S')’s and 1.4 x 10® r;6(3S)’s as described in the text while those in the 
e^e“ column are based on 10® T(35')’s assuming 250 fb“^ integrated luminosity. 


Parent Decay chain 


3® Si 


13.1% 


13.1% 3 


13.1% 3 


13.1% „3 


2®P27 (86.2) 1®D37 (96.5) 1®P27 (256.0) 

2®P27 (86.2) 1®P27 (104.4) ^ 1®P27 (248.4) 

2 ®P 27 (86.2) 

2 ®P 27 (86.2) 


1®P27 (104.4) 
0 . 02 % 3 


13.1% „3 


0.1 


13.1% 23 
13.1% o3 
12 . 6 °/ 
12 . 6 °/ 
12 . 6 °/ 
12 . 6 °/ 
12 . 6 '} 
5.9% 


2 ®P 27 (86.2) 

2®P27 (86.2) 

2 ®P 27 (86.2) 

2®Pi7 (99.3) 1®P27 (91.3) 

2®Pi7 (99.3) ^-1^” 1®P27 (91.3) 


3n_, 1103 31 i3p^^ (267.3) 

l®Pi7 (78.0) 1®P27 (239.1) 

l®Pi7 (78.0) ^ l®Pi7 (258.0) 


3 n... oi 411/0 ;^3p^^ (290.5) 


l®Pi7 (78.0) 
l®Pi7 (78.0) p+p- 

1.9% ., 3 , 


1)11? 135,^ (441.6) 2/12? p+p- 

l3s,^ (441.6) ^1^2? p+p- 

33.9% .,3„ i„oo ol 2.48% + _ 

->■ 1 Sl7 (423.0) —>■ p^p 

—)■ l®Si7 (441.6) —)■ p^p 

33.9% .,3„ /,r.o 2.48% + _ 

—)■ l®Si7 (423.0) —)■ p^p 

1-76% i3c- lool 11 2.48% + _ 

—>■ 1 Si7 (391.1) —>■ p^p 


2®Pi7 (99.3) 
2®Pi7 (99.3) 
2®Pi7 (99.3) 


0.80% .,3 


1®P27 (248.4) 
l®Pi7 (267.3) 
l®Pi7 (100.8) 1®P27 (239.1) 


l3s,^ (441.6) p+p- 

l3s,/y (423.0) ^122? p+p 

19.1% .3 o . 2.48% 


0.80% -3 


l®Pi7 (100.8) 225 l®Pi7 (258.0) 


33.9% ..3 


0.80% ..3 


rPi7 (100.8) 


47.1% ..3 


l®Po7 (290.5) 


1.76% ..3 


l®Si7 (441.6) 
l®Si7 (423.0) 


l®Si7 (391.1) 


p' p 


2.48% 


p^p 


2.48% 


PM 


3® So 


2®Po7 (122.0) 

2i25 2 ®Po 7 (122.0) 2 
2i25” 2 ®Po7 (122.0) 2 

1.8^-6 23^^^ (3Qg 

1.5x10“® ..3 


0.4% 


l®Pi7 (78.0) 2^” 1®P27 (239.1) l®Si7 (441.6) p+p 


4% 


4% 


19.1% . 3 , 


2.48% 


l®Pi7 (78.0) 225 l®Pi7 (258.0) 
l®Pi7 (78.0) "‘ll? l®Po7 (290.5) 


33.9% 3 


2.48% 


P P 


1.76% 3 


l®Si7 (423.0) 
l®Si7 (391.1) ^42|'” 


1.93% + _ 

2) —> p+p 


l®Si7 (840.0) p+p- 


Combined Events 


BR 

PP 

e"'"e 

6.8 X 10"® 

2100 

6800 

2.7 X 10"’’ 

84 

270 

1.6 X 10"® 

500 

1600 

2.0 X 10"® 

0.6 

2 

6.2 X 10"® 

19 

62 

5.4 X 10"® 

2 

5 

1.0 X 10"® 

0.3 

1 

2.5 X 10"® 

780 

2500 

1.5 X 10"® 

4650 

15,000 

7.6 X 10"® 

24 

76 

2.4 X 10"® 

740 

2400 

2.1 X lO"’’ 

65 

210 

1.8 X 10"® 

6 

18 

5.6 X lO"’’ 

170 

560 

4.8 X 10"® 

15 

48 

3.4 X 10"® 

5 

NA 

3.7 X lO"’’ 

52 

NA 
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TABLE XXVII: 4^S'i Decay chains, branching ratios and event estimates for LHCb Run II and Belle II. The number of events 
in the pp column are based on producing 2.3 x 10® T(45)’s as described in the text while those in the e+e“ column are based 
on 10^° T(4S)’s assuming 10 ab“^ integrated luminosity. The BR’s for the hadronic transitions were taken from the PDG |32| . 


Parent Decay chain 

Combined 

Events 


BR 

pp e''"e“ 


4®5i 


MM 


’ 3®Si 7 (171.6) p+p 


1.93% 


3®P27 (50.8) ’ 

3®P27 (50.8) 2®Si 7 (492.9) p+p' 

3®P27 (50.8) l®Si 7 (1013.8) p+p- 
3®P27 (50.8) 3®Si 7 (171.6) 2®Si7r+7r- p+p- 

3®P27 (50.8) 3®Si 7 (171.6) 


4.37% 3 


rSlTT+TT 


_|_ _ 2.48% _|. _ 


3®P27 (50.8) 2®Si 7 (492.9) I^Sitt+tt 

3®P27 (50.8) 2^Da 7 (72.7) 2^Pay (185.0) 

3^Pay (50.8) 2®P3 7 (72.7) 2®P27 (185.0) 

3®P27 (50.8) 


MM 

_ 2.48% + _ 

■ fl^fl 


2®Si7 (242.5) 

( 776 . 5 ) p+p- 


0^. 2®P3 7 (72.7) ^ 1®P27 (529.0) 


2.48% 


MM 


3®P27 (50.8) 


2^Da 7 (72.7) 1®P47 (96.7) 


100 % 


1^1)37 (200.9) 


91.0% .. 3 


RP27 (256.0) 


19 . 1 % i 3 


l®Si7 (441.6) 


2.48% 


MM” 


3^Pay (50.8) 
3®P27 (50.8) 


0.019% ..3 


0.13% „3 


1®P37 (350.0) 


2-^Da 7 (78.7) 


91.0% ..3 


56.2% 3 


1®P27 (256.0) 


2®Pi7 (191.6) 


19 . 1 % ,3 


19.9% „3 


RS17 (441.6) 


2®Si7 (229.6) 


2.48% 


MM 


1.93% 


MM 


7.2% 


2.18% 


MM 


3 ®Pi 7 (62.8) 3®Si 7 (159.8) 

3 ®Pi 7 (62.8) 2®Si 7 (481.4) p+p 


3 ®Pi 7 (62.8) l®Si 7 (1003.0) p+p 


2.48% 


3®Pi7 (62.8) 
3®Pi7 (62.8) 
3®Pi7 (62.8) 
3®Pi7 (62.8) ' 
3®Pi7 (62.8) ' 
3®Pi7 (62.8) 
3®Pi7 (62.8) ' 
3®Pi7 (62.8) ' 
3®Pi7 (62.8) 
3®Pi7 (62.8) ' 
3®Po7 (78.7) 
3®Po7 (78.7) 


7.2% 


7-2% o3 c. n rm ^3 s^n+n- 


2 . 6 % 


3®Si 7 (159.8) 
3®Si 7 (159.8) 
2® Si 7 (481.4) 
2®Pa 7 (66.8) 
2®Pa 7 (66.8) 
2®Pa 7 (66.8) 
2®P2 7 (66.8) 
2®P2 7 (66.8) 


2.82% 


235'i7r+7r 


_ 1.93% 


M+M~ 


17.85% 3 „ + _ 2.48% + _ 

- y 1 SlTT TT - > p p 


2-^Pay (178.6) 
2^Pay (178.6) 
2®Pi7 (191.6) 
2®Pi7 (191.6) 


2®Si7 (242.5) p+p- 


2®Pi 7 (74.7) ^ 2 ®Pi 7 (183.9) 


0.31% 3 


2®Pi 7 (74.7) 


3®Si 7 (144.0) 


l®Pi7 (541.2) 
2®Pi7 (183.9) 
2®Po7 (206.4) 


(776.5) 

3n''”2®Si7 (229.6) 

l3s,^ (764.3) p+p- 

435,^ (423.0) p+p- 


1.93% 


MM 


435,^ (743.4) p+p- 


2.18% 


MM 


^3 2 ^Si 7 (466.2) 

0.01% -3 


1.93% 


MM 


3 ®Po 7 (78.7) "224;” l®Si 7 (1003.0) 41+ 


3 ®Po 7 (78.7) 3®Si 7 (144.0) 

3 ®Po 7 (78.7) 3®Si 7 (144.0) I^Sitt+tt 


o 3 C + - 1 - 93 % + - 

2 OlTT TT - > M M 

_ 2.48% 


3 ®Po 7 (78.7) 2®Si 7 (466.2) l^Siir+rr- p+p 


2 . 2 ^- 33p^^ (78,7) o.(^% 

3®Po 7 (78.7) l®Pi7 (340.2) 


3p.., 1-78 '7^ 2 ®Pi7 (58.8) 115 2®Pi7 (183.9) 

1.6 


19 - 9 % o 3 ri /OOA 1 - 93 % + - 

—>• 2^5i7 (229.6) —> /i^/i 


1®P27 (239.1) 


2 .2x10" 


3®Po7 (78.7) 
3®Po7 (78.7) 


23 D_, 1 . 7 Q -71 o.oo9_i/„ ^ 3 ^^^ (340.2) l3pi^ (258.0) 


l®Si7 (441.6) p+p- 

4®Si7 (423.0) p+p- 


l®Pi7 (340.2) 


l®Po7 (290.5) 


i^^“4®Si7 (391.1) 


I X 10' 


3.1 X 10“® 7 

1.3 X 10“® 3 

1.0 X 10“® 2 

7.7x10“^° 0.2 

1.5 X 10“® 0.3 

2.9 X 10“® 0.7 

3.0 X 10“^° 0.07 

2.6 X 10“^° 0.06 

1.1 X 10“^° 


6.5 X 10“i^ 
3.0 X 10““ 
1.0 X 10“^° 

6.1 X 10“® 

1.9 X 10“® 
1.0 X 10“® 

1.5 X 10“® 
3.0 X 10“® 

4.4 X 10“® 

1.2 X 10“^° 
1.0 X 10“^° 
7.7 X 10“i° 

4.6 X 10“^° 

3.6 X 10“^° 

9.9 X 10“^^ 

9.5 X 10“^^ 

1.5 X 10“® 

3.3 X 10“^° 

5.5 X 10““ 

3.7 X 10““ 

7.4 X 10““ 

7.5 X 10““ 

6.5 X 10““ 

1.5 X 10““ 

4.7 X 10““ 
4.1 X 10““ 


0.03 

0.01 

0.02 

15 

4 

2 

0.3 

0.7 

1 

0.03 

0.02 

0.2 

0.1 

0.08 

0.02 

0.3 

0.08 

0.01 

0.02 

0.02 


)® 
310 
130 
100 
8 
15 

29 
3 

3 
1 

0.6 

0.3 

1 

610 

190 

100 

15 

30 
44 

1 

1 

8 

5 

4 
1 

0.1 

15 

3 

0.6 

0.4 

1 

1 

0.06 
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TABLE XXVIII: The 2P Decay chains, branching ratios and event estimates for LHCb Run II. These are based on producing 
1.0 X 10^ Xi)2(2P)’s, 9.1 X 10® Xf>i(2P)’s, 4.7 x 10® X6o(2P)’s and 7.4 x 10® hi,(2P)’s as described in the text. 


Parent Decay chain 


2®P2 


2® Pi 


2®Po 

2® Pi 


1®PI37 (96.5) 


1®P27 (256.0) 


(441.6) 

l3Six (441.6) ^312? mV- 


1®P27 (104.4) ^ 1®P27 (248.4) 

1345^7 (104.4) l®Pi7 (267.3) l®Si7 (423.0) m+P 


2.48% 


1®P27 (91.3) ^ 1®P27 (248.4) 
1®P27 (91.3) l®Pi7 (267.3) 


19.1% 3 


■\' 3 , n f i r ^\ 2.48% _|_ 

1 5i7 (441.6) — 
1^5i7 (423.0) 


2.48% 


0.80% 


l®Pi7 (100.8) 1®P27 (239.1) ^27” l®Si7 (441.6) 


li+p“ 


0.80% .,3 


l®Pi7 (100.8) 225 l®Pi7 (258.0) 


33.9% ..3 


l®Si7 (423.0) 


0.80% ..3 


l®Pi7 (100.8) ^-27° l®Po7 (290.5) 


1.76% 3 


1®5'i7 (391.1) 


0.4% 3 


l®Pi7 (78.0) 22| i3p2.y (239.1) l^S'n (441.6) 


19.1% 3j 


2.48% 


0.4% 3 


l®Pi7 (78.0) 225 l3pi7 (258.0) 


33.9% ..3 


l®Si7 (423.0) ^7l2|'» 


2 . 0 % 


48% 


22 % 


1®P27 (111-4) ^-27° l®Pi7 (262.5) ^ l®So7 (488.3) 
2®So 7 (257.7) 


49% 


100 % 


99 


100% 


1^5o 7 (825.8) 


100% 


gg 
' gg 


Combined BR Events 


5.2 X 10"® 

517 

2.1 X 10"® 

21 

1.2 X 10"® 

124 

2.0 X 10"® 

180 

1.2 X 10-^ 

1100 

6.1 X lO"’’ 

6 

1.9 X 10"® 

174 

1.6 X 10"® 

15 

3.0 X 10"’’ 

1 

9.4 X 10"® 

44 

9.0 X 10"® 

6.6 X 10‘ 

0.48 

3.6 X 10' 

0.22 

1.6 X 10' 
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TABLE XXIX: 3®Pj Decay chains, branching ratios and event estimates for LHCb Run II. These are based on producing 
9.9 X 10® Xi) 2 ( 3 P)’s, 7.9 x 10® Xf>i(3P)’s, 3.5 x 10® X6o(3P)’s and 7.3 x 10® /ii,(3P)’s as described in the text. 


Parent Decay chain 


3^P2 


3.8 

1.8 
1 . 1 % 


' 3^Si 7 (171.6) 


2^Si 7 (492.9) 
l^S"! 7 (1013.8) 


1.93% 




2.48% 


0.28% .3, 


2.48% 


1.93% 




MM 

l^PaTTTT l^Sn (441.6) —> M ' P 

“ill? l^Pirrir l®Si7 (423.0) mV" 

o^o ^3jj^ ^ (72.7) 2 ®P 27 (185.0) 

o^o 2375 ^ ^ (72.7) 2®P27 (185.0) l®S'i 7 (776.5) —> p ' ^ 

o^o ^ (72 7 ) i3p^^ (529.0) l®S'i 7 (441.6) /i+m" 

o^o ^3jj^ ^ (72 7 ) i3p^^ (gg 7 ) i_£^« (200.9) 


2®S'i7 (242.5) —> n 

2.48% 


19.1% 3 


0.019% .. 3 


13^37 (350.0) 1®P27 (256.0) 


19.1% 3 


1®P27 (256.0) 
l®S'i 7 (441.6) p+M 




1®5i7 (441.6) 


0.13% „3 


17% 


2®D2 7 (78.7) ^ 2®P27 (178.6) 


10 . 6 % 3 


0.13% 3 


17% 


7.0% 


0.13% „3 


56.2% 3 


2®Pi7 (191.6) 2^Siy (229.6) 
2®Pi 7 (191.6) l®S'i7 (764.3) 


0.13% „3 


0.13% 


2^S'i7 (242.5) 

2®P2 7 (78.7) 2®P27 (178.6) ^ l^Sij (776.5) 

2®D2 7 (78.7) 

2 ®D 2 7 (78.7) 

2 ®D 2 7 (78.7) 

2 ®P 2 7 (78.7) 


1.93% 


PM 


2.48% 


PM 


56.2% 3 


1.93% 


P P 


12 % 


3®Pi 


7.2% 


1 . 1 % 


3®S'! 7 (159.8) 

2®S'! 7 (481.4) /r+^i 


2.18% 

1.93% 


l^Pi 7 (541.2) 
13P37 (93.7) 

p+p- 


33.9% ..3 


l®S'i 7 (423.0) 
1®P27 (189.2) 


2.48% + _ 

—> P^M 
2.48% + _ 

—> P^M 

l®Pi7 (267.3) 

l®S'i7 (423.0) p+M" 


l^S"! 7 (1003.0) p+M 


0 . 75 % ,3, 


l®P27r7r l®Si 7 (441.6) 


2.48% 


P P 


0.48% 3, 


l®Si 7 (423.0) 


2.48% 


0.40% 3 


5.3x10“4% , _ 

-> P^M 


0.94% 3 


2®Pi 7 (74.7) 

2 ^D 2 7 ( 66 . 8 ) ^ 2 ®P 27 (178.6) 2 =^S'i 7 (242.5) /i+m 


10 . 6 % „3 ^ 


1.93% 




0.94% 3 


2®P2 7 ( 66 . 8 ) ^ 2®P27 (178.6) l^S'o (776.5) /r+^i 


7.0% 


2.48% 


0.94% ^3 


2'-^D2 7 ( 66 . 8 ) 2®Pi7 (191.6) 


19.9% „3 


2®S'i7 (229.6) 


1.93% 


+ /,- 


MM 


0.94% ^3, 


2®P2 7 ( 66 . 8 ) 2®Pi 7 (191.6) l®S'i 7 (764.3) fi+fi 


9.2% 


2.48% 


0.94% ^3 


23^2 7 ( 66 . 8 ) ^ 1®P27 (522.8) 


19 . 1 % .3 


l='S'i 7 (441.6) fi+fi 




0.94% ^3 


2®P2 7 ( 66 . 8 ) ^ l®Pi 7 (541.2) (423.0) /i+fi 


33.9% 3j 


2.48% 




O^o ^ 3 p^ ^ 6^0 ^ 3 p^^ 8^0 ^ 3 ^^^ (189.2) 


0.40% „3 


17% 


2®Pi 7 (74.7) ^ 2 ®Pi 7 (183.9) 


19.9% „3 


2^S'i7 (229.6) 


33 . 9 % .,3„ /.oo 2.48% j_ _ 

—s- l"^S^i 7 (423.0) —s- /i+ii 
1.93% 




0.40% „3 


2®Pi 7 (74.7) ^ 2®Pi 7 (183.9) l^S'o (764.3) 
2®Pi 7 (74.7) 


9.2% . 3 I 


0.40% „3 


2.48% 


MM 

p+p- 


28.1% „3 


2®Po 7 (206.4) ^214” 2®S'i7 (207.1) 


1.93% 


0.40% „3 


2®Pi 7 (74.7) 


28.1% „3 


2®Po 7 (206.4) l®S'i 7 (743.1) 


0.40% „3 


2®Pi 7 (74.7) 1 ®P 27 (90.5) 


2.48% 

74.7% 


MM 

■ fi+fj,- 


3i 7._. inn 71 i4.b/o (184.3) ^±14 i3p^.y (267.3 


3"Po 

3iPi 


0.054% 3 


l®P27r7r 


19.1% .,3 


1 ®S '17 ( 423 . 0 ) p+^- 


1®5'i7 (441.6) 


2.48% 


P P 


4.3% .. 1 „ nnoi nl 1007« 

—l^S'o 7 (1081.0) —> gg 

1 . 7 % .,1_ 49%.,1„ /.oo nl 100% 

—» VPi-K-K —» 1^5^0 7 (488.3) —» gg 


Combined BR Events 


8.3 X 10"^ 

8,240 

3.5 X 10"^ 

3,460 

2.7 X 10"^ 

2,710 

1.3 X 10"® 

132 

1.8 X 10"® 

173 

8.1 X 10"® 

80 

6.9 X 10"® 

68 

2.9 X 10"® 

29 

1.8 X 10"® 

18 

8.2 X 10"'^ 

8 

4.5 X 10"'^ 

4 

3.8 X 10"'^ 

4 

2.8 X 10"® 

28 

1.7 X 10"® 

17 

1.3 X 10"® 

13 

4.8 X 10"'^ 

5 

1.6 X 10"® 

12,360 

5.0 X 10"^ 

3,950 

2.7 X lO""' 

2,150 

3.5 X 10"® 

278 

4.0 X 10"® 

317 

2.1 X 10"'^ 

2 

3.3 X 10"® 

26 

2.8 X 10"® 

22 

2.0 X 10"® 

158 

1.2 X 10"® 

95 

8.9 X 10"'^ 

7 

1.1 X 10"® 

87 

3.5 X 10"® 

28 

2.6 X 10"® 

21 

1.6 X 10"® 

12 

1.0 X 10"® 

8 

2.5 X 10"'^ 

2 

1.4 X 10"'^ 

1 

2.6 X 10"® 

9 

4.3 X 10"® 

3.1 X 10 ' 

8.3 X 10"® 

6.0 X 10' 
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TABLE XXX: 4^Pj and 5®Pj Decay chains, branching ratios and event estimates for LHCb Run II. These are based on 
producing 9.7 x 10® Xf.2(4P)’s, 7.4 x 10® X6i(4P)’s, 3.1 x 10® X6o(4P)’s, 8.2 x 10® xm(5P)’s, 5.7 x 10® Xi>i(5P)’s and 2.2 x 10® 
Xf>o(5P)’s as described in the text. 


Parent Decay chain 






5®Pi 


5®Po 


0 . 012 % 3 1300 01 2 - 18 % + - 

—> S bi 7 (433.9) —> fd 


2^Si 7 (747.2) 


1.93% 


M M 


o.u A XU, ■ '/O w 3 ^ ^ - X 2.48% _|_ 

—^ 1 5i 7 (1255.1) —>• fi^fi 

335^ ^ (424.3) ^215? 




MM 


6 . 0 x 10 "^% ^ 3 ^ + _ 

—^ 2 bi 7 (737.9) —> ^ (I 


3.3x10" 


0-3 
0-3 
6.1 x10- ^ 

n-4 
> 

2.5x10-^ 


l^S"! 7 (1003.0) ^+^- 

^■4x10-’'” 3353 7 (411.8) 

l^S-i 7 (1234.8) fi+fi- 

ixio-% 335^ ^ (g4g g) 2^% 


3.4x10 


— 3q; 


’ 2®S'! 7 (953.7) fi+fi- 

' l^S"! 7 (1451.3) fi+fi- 

' 2^Si 7 (946.4) 

' l^Si 7 (1444.4) n+iJL- 

' 3®S'! 7 (629.9) /i+^- 

’ 2^Si 7 (937.3) /i+p- 

’ l^S-i 7 (1435.7) ^l+^JL- 


Combined BR Events 


2.6 X 10"® 

25 

2.6 X 10"'^ 

3 

1.3 X 10"® 

12 

2.7 X 10"® 

20 

1.1 X 10"'^ 

1 

8.2 X 10"'^ 

6 

1.4 X 10"® 

4 

1.5 X 10"'^ 

0.5 

2.1 X 10"® 

0.2 

4.8 X 10"'^ 

4 

8.4 X 10"'^ 

7 

2.5 X 10"'^ 

1 

3.5 X 10"'^ 

2 

6.5 X 10"® 

0.1 

1.3 X 10"'^ 

0.3 

1.0 X 10"'^ 

0.2 


TABLE XXXI: The IP Decay chains, branching ratios and event estimates for LHCb Run II and Belle 11. For the pp column 
these are based on producing 2.0 x 10'* T3(lP)’s and T3(lP)’s, 1.6 x 10* Ti(lP)’s and 1.9 x 10* i?62(lP)’s as described in 
the text. For the e''"e“ column these are based on 1.3 x 10® Ti(lP)’s produced assuming cr = 13 pb and 100 fb“* integrated 
luminosity. 


Parent Decay chain 


1®P3 

1®P2 


l®Pi 


1*P2 


rP27 (256.0) 

2.48% + 

—^ MM 

19.1% .,3 


0.74% - 3 ^ + 

- > 1 OlTT TT 

22 % 


435,^ (441.6) 


1®P27 (248.4) 43^^^ (444 g) 

l®Pi7 (267.3) l®S'i7 (423.0) p+p- 


0 . 66 % , 3 „ + _ 2.48% + _ 

- > 1 OlTT^TT - > P p 

l^Pa7 (239.1) 


28% 


^OQO (441.6) 


l®Pi7 (258.0) 1®5i7 (423.0) ^-2^” p+p 


2.48% 

2.48% 


MM 


47.1% - 3, 


l®Po7 (290.5) l®S'i7 (391.1) p+p 


2.48% 




0.40% 


2.48% 


l^S'lTT^TT - > /i ' fl 


3.93 x10- -^% + _ 

91:7 4 *Pi7 (262.5) ^ l*S'o7 (488.3) gg 

1 . 3 % , 1 „ + _ 100 % 

—1 Son- TT —)■ gg 


Combined Events 


BR 

PP 

e''"e 

4.3 X 10"® 

87 

NA 

1.8 X 10"* 

4 

NA 

1.0 X 10"® 

20 

NA 

6.3 X 10"® 

126 

NA 

1.6 X 10"* 

3 

NA 

7.6 X 10"® 

1 

99 

2.4 X 10"® 

38 

3120 

2.0 X 10"* 

3 

260 

9.9 X 10"® 

2 

129 

3.9 X 10"® 

1 

50 

0.45 

8550 

NA 

1.3 X 10"^ 

247 

NA 
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TABLE XXXII: The 2D Decay chains, branching ratios and event estimates for LHCb Run II. These are based on producing 
3.4 X 10'* T3(2 D)’s, 3.1 X 10* T2 (2D)’s, 2.4 x 10* Ti(2D)’s and 2.9 x 10* r/d2(2D)’s as described in the text. 


Parent Decay chain 


2^D3 


2* Da 


2^F2J (185.0) 


2^517 (242.5) 


65.1% „3 


2^P27 (185.0) l^S'iy (776.5) 


2.48% 




10.3% 3 


1=*P27 (529.0) 1^517 (441.6) m+m 


2.48% 


6.7% 


17% 


1=^P47 (96.7) 1®D37 (200.9) 

10.6% 3 j 


91.0% 3 


1*P27 (256.0) 


19.1% 3 


1 3 r-r / .1 ^ 1 /-‘\ 2.48% _|_ 

I Si7 (441.6) —>■ 


1.93% 


2 *Di 


2*P27 (178.6) 2*5'i7 (242.5) —> /r > 

^ 2®P27 (178.6) l®Si7 (776.5) 

2®Pi7 (191.6) 2^Si7 (229.6) p+M" 

5^0 23p^^ (191.6) i35'i7 (764.3) 

^ l®Pi7 (541.2) 135 'i 7 (423.0) /i+^i" 

1®P37 (93.7) 1^D27 (189.2) i3p^.y (267.3) 

5.28x10“3% 


33.9% t3„ /.ioo nl 2-48% + _ 

->■ 1 Pl7 (423.0) ->■ 


fX'fl 


17% 


2^517 (229.6)'^’'VV- 
*_:^»235l7(207.1)^^’'W 
o_:55 13537 (743.1) 


2* Da 


67.1% 


2®Pi7 (183.9) 

2^Po7 (206.4) 

2^Po7 (206.4) 

1®P27 (90.5) 

2*Pi 7 (188.3) ^ 2*Po7 (257.7) 


1^D27 (184.3) 


100 % 


l^Pi7 (267.3) 
99 


33.9%,3„ /.r,no3 2.48% J. _ 

—> l*Pi7 (423.0) —> /i+fj, 


67.1% „i 


2*Pi 7 (188.3) ^ l*Po7 (825.8) 


100 % 


99 


12 % 


l*Pi7 (536.5) 1*507 (488.3) 


100 % 


99 


Combined BR Events 


1.3 X 10"^ 

45 

1.1 X 10"^ 

38 

4.9 X 10"* 

16 

2.9 X 10"* 

10 

3.5 X 10"* 

11 

3.0 X 10"* 

9 

2.1 X 10"® 

66 

1.3 X 10"^ 

39 

1.0 X 10"® 

31 

3.7 X 10"* 

11 

5.3 X 10"® 

1 

6.5 X 10"* 

16 

2.5 X 10"* 

6 

6.3 X 10"® 

2 

3.6 X 10"® 

1 

0.32 

9280 

0.15 

4350 

0.059 

1710 


TABLE XXXIII: The IP Decay chains, branching ratios and event estimates for LHCb Run 11. These are based on producing 
166 X54(lP)’s, 163 Xf,3(lP)’s, 147 X52(lP)’s and 158 /if,3(lP)’s as described in the text. 


Parent Decay chain 

1®P4 l^Dax (184.9) l^Pay (256.0) l^Sij (441.6) pi+fi- 

l3p3 13^3^ (189 2) (267.3) l^Siy (423.0) ^i+m" 

l^Pa *■24'° l^Dix (194.1) 24 l^Pi7 (258.0) *24° l®5i7 (423.0) *-24 ^+/i- 
1*P3 (4° 1*D27 (189.2) '^24° l*Pi7 (262.5) ^4 l*5o7 (488.3) *-4° gg 


Combined BR Events 


4.3 X 10"® 

0.7 

5.6 X 10"® 

1 

1.9 X 10"® 

0.3 

0.45 

71 
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TABLE XXXIV: The 2^_Di Decay chains, branching ratios and event estimates for Belle II. The event numbers are based on 
2 X 10® T3(2D)’s produced assuming a = 18 pb and 100 fb“^ integrated luminosity. 


Parent Decay chain 


Combined BR Events 




5.28x10“-^% + _ 

-1 M ft 


1 % 


17.2% „3 


2.4% 


28.1% „3 


28.1% 3 


2®P27 (170.9) 
2®P27 (170.9) 

” 1^P27 (515.4) 
° 23Pi7 (183.9) 
23Pi7 (183.9) 
l®Pi7 (533.8) 
23Po7 (206.4) 
2 ®Po7 (206.4) 


2®5i7 (242.5) 

135,^ (776.5) mV 

135^7 (441.6) "-^’'W 
l^^»2®5l7(229.6)^^^W 

9.2% 


5.3 X 10"® 
2.0 X 10"® 
1.7 X 10"® 

2.4 X 10"® 

6.6 X lO""* 


33.9% 
->• 

4.6% 


0.9% 


7.7% 


l®Po7 565.4) 
1®P27 (90.5) 
1®P27 (90.5) 
1®P27 (90.5) 


1.76% 


1®D27 (184.3) 
1®D27 (184.3) 
1®D37 (176.5) 


1®P27 (90.5) 1®Di7 (194.1) - 


1®P27 (90.5) 
1®P27 (90.5) 


1®Di7 (194.1) 
1®Di7 (194.1) 


l®Pi7 (267.3) 
1®P27 (248.4) 
1®P27 (256.0) 
’ 1®P27 (239.1) 
l®Pi7 (258.0) 
l®Po7 (290.5) 


31:5? 13537 (423.0) 

1 ) 51 ? 135 ^^ ( 441 . 6 ) 3 : 12 ? ^+^- 

151? 135^7 (441.6) 

115 ? 135 ^^ ( 441 . 6 ) 315 ? ^+^- 

315? 135,7 (423.0) ^A’'VV- 
ll5?l35^.,(391.1)^-55?pV 


3.7 X 10"® 
6.1 X 10"® 


106 

40 

34 

5 

1320 


^ i35i 7 (764.3) /i+^- 

3.9 X 10-^ 

780 

' l®^^ (423.0) /i+^- 

2.0 X IQ-^ 

400 

■ 23517 ( 207 . 1 ) 

2.5 X 10-^ 

500 

■ l®^^ (743.1) 

6.3 X 10"® 

126 

l®^^ (391.1) '*-25?/i+^“ 

3.4 X 10"® 

68 


74 

12 

1 

5 

162 

14 


TABLE XXXV: The 3^Di Decay chains and combined branching ratios. These are small enough that we do not expect enough 
events at Belle II for this state to be observed. 


Parent Decay chain 


33 Di 


2.30x10“®% j. _ 
-> 

5.6x10“ ' 


33P27 (168.6) 


3.8% „3„ 2.18% + _ 

—>• 3 5i7 (171.6) —>• 


5.6x10 

9.2x10^ 

9.2x10' 


3^P27 (168.6) 
33Pi7 (180.4) 


15? 23517 (492.9) 


3d_, 7180 41 335 ,^ (159.8) 


2.18% 


M M 


1.35x10“"'% o3 


33Pi7 (180.4) 
33Pi7 (180.4) 

-2 (y 

33Po7 (196.2) 
33Po7 (196.2) 


33Po7 (196.2) 


5 : 5 ? 23517 ( 481 . 4 ) 

l^Siy (1003.0) ^-25? ^+^- 
» 255 ? 335,^ ( 144 . 0 ) 3 : 55 ? ^ V" 

23517 (466.2) ^755? ^+^- 
0.01% ,3 


2.48% 




i3S'i7 (988.5) —> M' fJ- 
1.93% 


' 23P27 (420.4) 23517 (242.5) 


2 x 10 - 5 % 3 


23P27 (196.2) 1^” 135i7 (776.5) 


2.48% 




9.3x107**% ^3 


19-9% ^3 


2*^517 (229.6) 


9.3x10““*% ^3 


2^Pi7 (433.8) 


1.93% 


MM" 


23Pi7 (433.8) i3S'i7 (764.3) 


2.48% 


MM 


1.3x10““*% ., 3 


33.9% -3 


rSij (423.0) 


2.7x10“'^% 3 


rPij (533.8) 


2.48% 


MM 


23Po7 (454.9) 23517 (207.1) 


1.93% 


/7 


2.7x10“®% 3, 


23Po7 (454.9) l^Siy (743.1) /r+ii 


2.48% 


5.7x10“'*% 


l3Po7 (806.1) 
23p27 (75.7) 


1.76% 3 


i3S'i7 (391.1) 


2.48% 


‘ 2^Dij (172.6) 


PM 

5.28x10" 


Combined BR 

2.3 X 10"3 
4.6 X 10"® 

1.9 X 10"® 

1.4 X 10"'^ 

4.6 X 10"3 

2.5 X 10"3 

9.1 X 10"® 

2.0 X 10"® 

3.3 X 10"^® 

4.1 X 10"^® 

3.5 X 10"^® 

3.6 X 10"® 

2.1 X 10"® 

1.1 X 10"® 

2.4 X 10"® 

6.0 X 10"® 

2.5 X 10"® 


/J.+/I- 2.7 X 10” 
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TABLE XXXVI: The 4^_Di Decay chains and combined branching ratios. These are small enough that we do not expect enough 
events at Belle II for this state to be observed. 


Parent Decay chain 


Combined BR 


. 3 _, 3.04x10“'’% j. _ 

4x10“^% 3 


4x10“^% 

n-5 
> 

1.8x10-3% „3 


3^P27 (392.7) 3®S'i7 (171.6) 

3®P27 (392.7) 

^P 27 (392.7) 

33Pi7 (404.2) 


PM 

2®S'i 7 (492.9) p+M" 
3^P27 (392.7) 1^5r7 (1003.0) 


n~3Q 


3®5'i7 (159.8) — 
235'i 7 (481.4) —> 
i.8x_m^3% (404.2) l^S'o (1003.0) 


3'^Pi7 (404.2) 

3®Pi7 (404.2) 

5.6x_^- 3% 33p^^ (4^9 g) 0^% 33 ^^.^ 2 ^% 


PM 

1.93% + _ 


5.6 xj^ 3% (419.6) o.^% 23^^.^ (466.2) 

5.6xj^3% ^ 3^^^ (419.6) i3s^y (988.5) /i+m" 

2.ixj^3% 23p3^ (639.1) 2^S'i7 (242.5) /i+m" 

1x^3% ^3 ^^^ (652.3) 2^517 (229.6) ^+ai- 

ixi^"% (652.3) l®S'i7 (764.3) p+M" 
3.6xj^‘‘% 23p^^ (672.9) 235 'i 7 (207.1) /i+n” 

3.6x_^-^% 23p^^ (072.9) 13537 (743.1) /r+p- 

2.6x_^-^% ^3p^^ (986.0) 13537 (423.0) /iV 

i.ox_io^-3% ^3p^^ ^016.7) 1^517 (391.1) V 


0 

CO 

X 

10"® 

3.3 

X 

10-10 

1.4 

X 

10-10 

1.1 

X 

10-10 

2.8 

X 

10"3 

9.0 

X 

10-9 

4.9 

X 

10-9 

bo 

X 

10-9 

00 

CO 

X 

10-10 

1.4 

X 

10-10 

4.3 

X 

10"“ 

bo 

X 

10-9 

2.3 

X 

10-9 

3.2 

X 

10-9 

00 

0 

X 

10-10 

2.2 

X 

10"® 

4.4 

X 

10-9 



